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A  class  of  multiple  facility,  multiple  product,  production- 
inventory  scheduling  problems  is  considered  over  a  finite  planning 
horizon.   The  horizon  consists  of  discrete  production  periods  during 
each  of  which  at  most  one  product  can  be  assigned  to  a  facility. 
Product  demands  are  constant  over  a  period  but  not  necessarily  identi- 
cal in  all  periods.   Each  product  has  a  minimum  batch  size  which  can 
be  made,  and  only  integer  multiples  of  this  basic  batch  size  are  allowed. 
The  problem  objective  is  to  determine  an  assignment  of  products  to  facil- 
ities which  minimizes  the  sum  of  any  production,  inventory,  backordering 
and  facility  usage  charges  which  occur  over  the  horizon.   Four  differ- 
ent scenarios  for  this  class  of  problems  are  considered. 

In  a  straightforward  manner,  each  of  the  scenarios  which  are 
considered  is  formulated  as  an  integer  program.   It  is  then  shown  that 
each  of  the  problems  can  be  reformulated  as  a  network  flow  problem  which 
can  be  solved  using  any  of  the  very  efficient  algorithms  for  finding 
minimal-cost  flows  in  single  commodity  networks.   In  addition,  for  some 


important  special  cases,  single  pass  procedures  are  developed  which  are 
more  efficient  than  the  standard  flow  algorithms.   In  the  fourth  sce- 
nario the  horizon  length  is  also  considered  to  be  a  variable  in  the 
optimization.   In  this  scenario  a  specially  structured  integer  program- 
ming formulation  is  identified,  and  an  extremely  efficient  algorithm  is 
developed  for  determining  its  solution. 


CHAPTER  1 


INTRODUCTION 


An  Overview 


Inventory  management  and  production  control  represent  two  major 
components  of  the  manufacturing  process.   Traditionally,  research  in 
production-inventory  systems  has  been  divided  into  two  groups  based  on 
the  assumptions  regarding  the  interaction  of  these  components.   In  one, 
the  assumption  is  made  that  the  two  components  are  independent,  and 
optimization  can  be  achieved  by  optimizing  the  control  of  each  component 
separately.   In  the  other,  inventory  management  and  production  control 
are  not  assumed  to  be  independent,  but  rather  the  components  interact 
and  optimization  must  be  achieved  by  considering  the  two  components 
jointly. 

The  research  in  this  dissertation  falls  in  the  second  group  of 
thought  and  deals  more  specifically  with  the  analysis  of  a  multi- 
facility,  multi -product,  integrated  production-inventory  system. 
Realistic  examples  of  the  multiple  facility,  multiple  product  system 
include  textile  spinning  machines,  metal  stamping  operations,  steel 
rolling  facilities,  glass  molding  lines,  container  manufacturing  machines 
and  packaging  lines. 

Contributions  within  this  dissertation  include  original  formu- 
lations of  the  multi-facility,  multi-product  problem  as  integer  programs 


and  network  flow  problems  which  offer  considerable  insights  into  the 
operation  of  the  production-inventory  system.   These  insights  allow  us 
to  enrich  the  problem  generously  to  include  many  other  realistic  aspects 
of  the  production-inventory  system  heretofore  not  solvable  with  present 
techniques.   In  addition,  efficient  algorithms  are  developed  for  solving 
different  types  of  associated  problems  which  in  many  cases  are  single 
pass  procedures  which  can  be  easily  carried  out  by  hand. 

Literature  Review 

Problems  of  the  type  considered  in  this  dissertation  generally 
have  been  modeled  as  one-facility,  deterministic  problems  where  produc- 
tion rates  and  demand  rates  for  each  of  the  products  take  on  known  con- 
stant values.   Only  one  product  can  be  produced  on  the  facility  at  a 
time,  and  a  setup  is  necessary  each  time  a  different  product  is  to  be 
produced.   Stockouts  are  not  permitted,  and  backordering  to  meet  demands 
may  or  may  not  be  allowed.   Inventory  carrying  charges  for  each  of  the 
products  are  assessed  in  direct  proportion  to  the  amount  of  product  on 
hand.   The  problem  is  to  determine  a  schedule  for  allocating  the  produc- 
tion time  of  the  facility  to  each  of  the  products  so  that  the  costs  asso- 
ciated with  setting  up,  producing,  backordering  and  carrying  inventory 
are  minimized  over  a  given  planning  horizon. 

One  general  difference  in  models  dealing  with  this  problem  is 
in  the  definition  of  the  planning  horizon.   In  one  approach  the  length 
of  the  planning  horizon  is  assumed  to  be  infinite,  and  in  the  other  the 
length  of  the  horizon  is  assumed  to  be  fixed  and  finite.   In  the  finite 
horizon  case,  the  optimization  problem  involves  the  determination  of 


3  feasible  assignment  of  products  to  the  facilities  which  minimizes 
costs.   In  the  infinite  horizon  approach,  the  existence  of  a  finite, 
but  variable  length,  production  cycle  is  assumed  where  this  production 
cycle  repeats  itself  infinitely  many  times  throughout  the  horizon. 
Hence,  the  infinite  horizon  problem  involves  the  joint  determination 
of  the  optimal  production  cycle  and  the  optimal  assignment  of  products 
to  the  facility  during  this  production  cycle. 

A  well-known  approach  to  the  infinite  horizon  form  of  the  problem 
is  to  assume  that  each  product  is  set  up  and  produced  exactly  once  dur- 
ing each  production  cycle  where  the  amount  produced  for  each  product  at 
each  setup  is  proportional  to  the  length  of  the  production  cycle.   In 
this  case  the  object  of  the  optimization  reduces  simply  to  the  determi- 
nation of  the  optimal  length  of  the  production  cycle.   This  approach  is 
typically  referred  to  as  a  model  with  a  pure  rotation  cycle.   An  optimal 
macro-plan  for  this  model  has  been  determined  by  Magee  [12]  and  by 
Rodgers  [14]. 

The  more  general  approach  to  the  infinite  horizon  problem  is 
not  to  require  a  pure  rotation  cycle.   Rather  each  product  may  be  pro- 
duced on  an  individual  production  cycle  where  the  only  restriction  is 
that  feasibility  must  be  maintained  with  respect  to  the  facility.  At 
present,  there  is  no  general  method  for  determining  an  optimal  macro- 
plan  for  this  model.   However,  upper  and  lower  bounds  can  be  determined 
on  the  cost  of  the  optimal  solution.   A  lower  bound  is  determined  by 
solving  the  one-facility,  one-product  problem  with  a  pure  rotation 
cycle  for  each  of  the  products  and  summing  their  respective  costs. 
If,  in  fact,  the  optimal  production  quantities  determined  by  solving 


each  of  these  problems  can  be  scheduled  on  the  facility  without  causing 
scheduling  conflicts,  then  this  schedule  is  an  optimal  macro-plan  for 
the  multiple  cycles  problem.   In  like  manner,  an  upper  bound  on  the 
optimal  cost  is  determined  by  solving  the  one-facility,  multiple  product 
problem  with  a  pure  rotation  C3'cle  and  noting  the  associated  cost. 
Researchers  who  have  investigated  this  problem  include  Bomberger  [1] , 
Elmaghraby  [3,41,    Haessler  [9],  Hodgson  [10],  Madigan  [11],  Magee  [12], 
Maxwell  [13],  and  Stankard  and  Gupta  [15]. 

Of  primary  interest  in  this  dissertation  is  the  finite  horizon 
approach  to  the  problem.   In  this  approach  the  horizon  is  assumed  to 
consist  of  discrete  production  periods  (shifts,  for  example)  during 
which  at  most  one  product  can  be  assigned  to  the  facility.   During  any 
period  product  demands  are  assumed  to  be  constant  but  not  necessarily 
the  same  in  all  periods. 

The  production-inventory  scheduling  problem  was  initially 
modeled  as  a  finite  horizon  problem  by  Wagner  and  Whitin  [16].   They 
considered  a  one-facility  problem  for  the  case  where  backordering  is 
not  allowed  and  production  amounts  in  any  period  can  vary.   They 
developed  a  dynamic  programming  procedure  for  solving  the  problem  and 
were  able  to  realistically  enrich  the  problem  to  include  period  depend- 
ent product  demands,  inventory  costs  and  setup  costs. 

Zangwill  [17]  extended  the  results  of  Wagner  and  Whitin  to 
include  backordering  and  showed  that  the  problem  could  be  formulated 
and  solved  as  a  minimal-cost  flow  problem.   In  addition,  he  generalized 
the  cost  function  to  include  concave  inventory  and  backordering  costs. 


In  order  to  bring  the  finite  horizon  approach  to  the  problem 
more  closely  in  parallel  to  the  real  environment,  Elmaghraby  and  Mallik 
[2  ]  considered  the  problem  in  a  slightly  different  context.   They 
assumed  that  the  production  amounts  of  each  product  at  each  assignment 
of  the  product  to  the  facility   are   fixed  and  constant.   They  then 
developed  a  dynamic  programming  recursion  for  solving  the  problem  for 
the  case  where  demand  rates  and  production  costs  are  assumed  to  be  con- 
stant and  identical  in  all  periods.   For  the  batch  type  production  sys- 
tem where  product  inventories  are  updated  only  at  the  end  of  each  pro- 
duction period,  their  results  apply  to  a  special  case  of  the  problem 
considered  in  Chapter  3  of  this  dissertation. 

Practical  Motivations 

While  the  assumptions  in  the  previous  models  provide  a  basis  for 
analyzing  industrial  scheduling  problems,  it  is  frequently  true  that  the 
characteristics  of  the  real  system  are  much  more  complex  than  these 
models  allow.   For  example,  many  production  operations  consist  of  numer- 
ous facilities  (M) .   When  customer  demands  are  received  which  require 
processing  in  this  operation,  the  system  controller  must  decide,  in 
addition  to  what,  when  and  how  much  of  the  product  to  produce,  which 
facility  or  facilities  should  be  used.   Underlying  in  the  previous 
models  is  the  assumption  that  facility  selection  is  a  deterministic 
problem  which  is  not  part  of  the  decision  making,  and  hence  system 
optimization  can  be  achieved  by  solving  M  one-facility  problems. 
However,  this  is  not  always  the  case  in  the  real  environment  as  system 


controllers  frequently  interchange  products  between  facilities  in  order 
to  meet  demand  requirements  and  minimize  costs. 

Another  consideration  is  initial  product  inventories.   In  the 
infinite  horizon  approaches,  the  assumption  is  made  that  the  initial 
product  inventories  are  "in  phase  with"  the  optimal  macro-plan  of  allo- 
cating the  production  time  of  the  facility  to  each  of  the  products. 
In  the  industrial  environment  this  is  not  necessarily  true,  and  indus- 
trial schedulers  in  using  the  infinite  horizon  approach  may  still  be 
left  with  a  subhorizon  of  unknown  time  duration  during  which  the  facil- 
ity must  be  scheduled  "out  of  phase"  with  the  optimal  macro-plan. 

In  the  industrial  environment,  it  is  rare  that  decision  making 
is  based  on  a  single  plan  of  operation  as  many  of  the  previously 
described  models  assume.   For  example,  an  industrial  scheduler  fre- 
quently may  vary  production  run  lengths,  schedule  overtime,  schedule 
alternate  facilities,  purchase  product  from  outside  sources,  etc.,  in 
order  to  effectively  control  the  system. 

Another  consideration  which  is  frequently  applicable  to  many 
production  operations  is  a  fluctuating  demand  rate.   The  fluctuations 
of  the  demand  rate  can  arise  in  numerous  ways  at  least  some  of  which 
can  be  predicted  in  advance  by  good  forecasting.   Many  industries  have 
seasonal  demands  year-in-and-year-out  where  the  demand  might  be  repre- 
sentable  by  a  cyclic  function.   Some  industries  have  a  record  of  steady, 
continued  growth  which  can  be  represented  possibly  by  a  ramp  or  exponen- 
tial growth  function.   Other  industries  have  consistent  sales  through- 
out the  year  with  large  increases  around  the  holiday  season  (impulse 
function).   Temporary  plant  shutdowns  for  maintenance  purposes  or 


vacations  also  occur  in  which  case  the  demand  rate  for  the  supplier  to 
these  plants  may  suffer  a  sudden  drop  (negative  impulse  function).   New 
plants  are  being  bu .    and  new  uses  are  being  found  for  different 
products  either  of   r>ich  might  correspond  to  a  step  function  in  the 
demand  rate  forecast. 

As  these  examples  indicate,  the  real  environment  is  often  much 
more  complex  than  the  assumptions  of  existing  models,  and  the  develop- 
ment of  new  insights  is  necessary  if  optimal  system  control  is  to  be 
achieved.   This  dissertation  deals  specifically  with  some  of  the  prob- 
lems associated  with  the  operating  production-inventory  system.   The 
objective  of  the  research  is  to  provide  original  insights  into  the 
achievement  of  optimal  control  in  this  complex  system  as  well  as  to 
provide  efficient  procedures,  based  on  these  insights,  for  solving 
related  problems  in  the  real  environment. 

A  General  Problem  Description 

An  industrial  process  made  up  of  M  facilities  in  parallel  is 
considered  where  there  are  N  different  products  to  be  produced  over 
a  finite,  but  possibly  variable  length,  planning  horizon.   The  horizon 
is  made  up  of  production  periods  (shifts,  for  example)  which  have  equal 
time  duration.   During  any  period  at  most  one  product  can  be  scheduled 
on  each  facility.   When  any  product  i  is  assigned  to  a  facility  for  one 
period,  an  integer  multiple  of  p.  units  of  product  is  made.   Production 
completed  is  added  to  each  product's  inventory,  and  nonnegative  demands 
(d   units  for  product  i  in  period  k)  are  satisfied  from  this  inventory. 


8 

Backordering  may  or  may  not  be  allowed  (depending  on  the  application) , 
and  setups  are  assumed  to  occur  between  production  periods.   Initially, 
each  product  i  is  assumed  to  have  a  nonnegative  beginning  inventory 
level  denoted  by  I    (units  of  product  i) ,  and  when  applicable,  a  non- 
negative  initial  backorder  level  denoted  by  B    (units  of  product  i) . 
During  the  horizon  the  system  controller  must  schedule  the  products  to 
the  facilities  in  such  a  manner  as  to  meet  all  product  demands,  and  at 
the  same  time  build  up  each  product' s  ending  inventory  level  to  a 

desired  nonnegative  amount  denoted  by  I.„.   The  problem  objective  is  to 

iH 

determine  such  a  production  schedule  which  minimizes  the  sum  of  any  pro- 
duction, inventory,  backordering  and  facility  usage  charges  over  the 
horizon. 

As  an  example  of  this  problem,  consider  a  simplified  version  of 
a  two-product,  one-facility  problem  over  a  three-period  horizon  where 
at  most  p.  units  of  each  product  i  can  be  produced  by  the  facility  dur- 
ing any  period.   During  any  period,  then,  the  system  controller  has 
three  options.   He  can  schedule  either  product  one  or  product  two  to 

the  facility,  or  he  can  leave  the  facility  idle.   Hence  in  this  problem 

3 

there  are  3   or  27  possible  solutions  to  the  scheduling  problem.   Depend- 
ing on  the  imposed  constraints,  all,  some  or  none  of  these  schedules 
may  be  feasible  solutions  to  the  system  controller' s  problem.   What  the 
system  controller  must  attempt  to  do  is  determine  a  feasible  schedule 
which  minimizes  costs.   In  general  for  a  N-product,  M-facility,  H-period 

VTTJ 

problem  of  this  type,  there  will  be  (N+l)    possible  schedules  to  eval- 
uate.  Hence,  it  is  clear  that  something  significantly  better  than  total 
enumeration  is  needed  for  the  system  controller  to  be  able  to  solve 


realistic  kinds  of  problems.   One  of  the  goals  of  this  dissertation  is 
to  provide  an  efficient  procedure  for  solving  this  problem. 

Organization  of  the  Dissertation 

In  Chapter  2  a  fixed  horizon  problem  with  production  and  inven- 
tory costs  is  addressed.   A  mathematical  formulation  of  this  problem  is 
developed  which  allows  the  problem  to  be  solved  using  any  of  the  algo- 
rithms for  finding  minimal-cost  flows  in  single  commodity  networks. 
However,  for  an  important  special  case,  a  single  pass  algorithm  is 
given  which  is  even  more  efficient  than  the  standard  flow  algorithms. 

In  Chapter  3  a  fixed  horizon  problem  with  production,  inventory, 
and  backordering  charges  is  considered.   In  a  straightforward  manner, 
this  problem  is  formulated  as  a  linear,  mixed  integer  program.   It  is 
shown,  however,  that  the  problem  can  be  reformulated  as  an  all  integer 
program  which  can  be  solved  as  a  minimal -cost  flow  problem. 

In  Chapter  4  a  fixed  horizon  problem  is  considered  where  the 
system  controller  has  an  overtime  option  available  for  use  in  schedul- 
ing the  system.   This  problem  is  formulated  as  an  integer  program,  and 
a  solution  procedure  is  developed  for  solving  the  problem  which  involves 
solving  a  minimal-cost  flow  problem. 

In  Chapter  5  the  problem  of  bringing  initial  product  inventories 
'in  phase  with"  a  desired  infinite  horizon  cycle  plan  is  considered. 
The  problem  is  formulated  as  an  integer  program  which  has  special  struc- 
ture, and  a  very  efficient  algorithm  is  developed  for  solving  the  problem. 

A  summary  of  the  results  of  this  research  and  suggestions  for 
future  research  are  given  in  Chapter  6. 


CHAPTER  2 


A  NETWORK  APPROACH  TO  THE  PRODUCTION-INVENTORY 
SCHEDULING  PROBLEM  WITHOUT  BACKORDERING 


Introduction 

The  production-inventory  scheduling  problem  is  considered  over 

a  fixed  planning  horizon  made  up  of  H  production  periods.   Each  of  the 

M  facilities  is  considered  to  be  identical  in  that  product  i  can  be 

produced  by  any  facility  at  a  fixed  production  rate  of  p.  units  per 

period  and  a  cost  of  c.  dollars  per  period.   Demands  for  each  product  i 

occur  continuously  during  each  period  k  according  to  a  fixed  demand 

rate  of  d   units  per  period.   Inventory  carrying  charges  are  assessed 

on  product  inventories.   During  the  horizon  the  system  controller  must 

schedule  the  products  to  the  facilities  in  such  a  manner  as  to  meet  all 

demands,  without  backorders ,  and  at  the  same  time  build  up  each  product 

i's  inventory  level  to  a  desired  minimum  ending  amount  denoted  by  I.  . 

iH 

The  objective  is  to  determine  a  production  schedule  which  minimizes  the 
sum  of  production  and  inventory  carrying  charges  over  the  horizon. 

This  problem  is  formulated  as  a  network  flow  problem  for  which 
efficient  solution  procedures  are  known.   However,  for  an  important 
special  case  a  single  pass  algorithm  is  given  which  is  more  efficient 
than  the  standard  flow  algorithms.   In  addition,  the  general  model  can 
be  enriched  generously  to  include  numerous  other  realistic  aspects  of 

10 
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the  production  system  without  severely  complicating  the  solution  proce- 
dure.  These  generalizations  include  such  things  as  product  and  facil- 
ity dependent  production  costs,  facility  usage  costs,  period  dependent 
facility  availability,  period  dependent  inventory  carrying  costs,  inven- 
tory constraints,  facility  and  product  assignment  restrictions,  and 
horizon  dependent  facility  allocation  costs.   For  presentation  purposes, 
only  basic  assumptions  are  used  in  the  development  of  the  problem  with 
the  more  detailed  aspects  of  the  model  being  identified  in  the 
"Extensions  of  the  Model"  section.   Many  of  these  enrichments  will  no 
doubt  become  apparent  in  the  development. 

Problem  Development  and  Formulation 

For  the  problem  being  considered,  backordering  is  not  allowed 
and  product  inventory  changes  during  any  period  are  assumed  to  occur 
linearly  (i.e. ,  during  any  period  the  production  and  demand  rates  are 
constant).   Thus  stockouts  on  product  i  will  occur  during  period  k  if 
and  only  if  at  least  one  of  the  ending  period  inventories,  I.      or 
I   ,  is  negative.   Hence,  requiring  that  ending  period  inventories  be 
nonnegative  for  all  products  and  periods  is  a  necessary  and  sufficient 
condition  for  preventing  stockouts. 

Define  w   as  the  minimum  number  of  times  (integer)  product  i 
must  be  produced  during  the  first  k  periods  in  order  to  have  a  non- 
negative  inventory  level  at  the  end  of  period  k.   Then 


k 

ij     iO   '  i  J     k  =  1,2,... ,H-1 


{o,  «5   dj,  -i,„>/p>      '  ~  1-2---"s 


d=i 
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where  (a)  denotes  the  smallest  integer  greater  than  or  equal  to  a. 

In  addition  to  being  nonnegative,  product  i's  inventory  level  at  the 

end  of  the  horizon  must  be  greater  than  or  equal  to  I   .   Hence,  define 

iH 


H 
{0,   <(I    +   E   d    -  Ii0)/Pi>3     1-  1.2,. ..,N, 
j=l 


and  let  W  be  an  N  by  H  matrix  with  components  [w   ]. 

Define  x.,  as  the  number  of  facilities  scheduled  to  produce 
lk 

product  i  during  period  k  where  X  is  defined  as  an  N  by  H  matrix  with 
component  [x   ].   The  inventory  level  for  product  i  at  the  end  of  per- 
iod k  can  then  be  written  as 


I.,=P.   £  x.  .   +    I.n   -     Ed..,  k  =  1,2,...  ,H. 

x  j=i  XJ  l0   j=1  XJ  • 


Figure  2.1  gives  an  example  illustration  of  the  change  in  inventory  of 

product  i  during  a  period  in  which  the  amount  produced  on  product  i 

exceeds  the  amount  demanded.   Obviously,  the  linear  change  in  inventory 

of  product  i  during  period  k  is  a  function  of  the  number  of  facilities 

scheduled  to  produce  product  i,  and  the  slope  of  the  line  can  take  on 

M+l  possible  values  (i.e.,  x   can  be  0,  1,2,...,M-1  or  M)  in  period  k. 

lk 

The  average  amount  of  inventory  on  hand  for  product  i  during  period  k 
in  any  of  these  cases  is  given  by 


k-1  k-1 

(-){l.    .     _  +  I..]=  (-)fp.x.,  +  2E  x.  .  +  21.  _  -d      -2   E  d.  .} 
2        i,k-l        ikJ         2        i    lk       .    „    ij  lO        lk       ,   „   ijJ 

k  =    1,2,. . . ,H 

where  x.^  and  d.„  are  defined  to  be  zero. 
lO      lO 
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i,k-l 


TIME  (IN  PERIODS) 


Figure  2.1   Product  Inventory  Change  during 
a  Production  Period. 
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Denoting  by  0.  the  inventory  carrying  charge  per  unit  of  product 
i  carried  in  inventory  for  one  period,  where  0.  >  0,  the  total  inventory 
cost  over  the  horizon  is  given  by 

N   H  1  k-1       N   H  k-1 

E   E  (-)0. p. {x   +  2  E  x.  .}+  E   E  (-)0.{2I.  -d.  -2  E  d.  .}. 
i=i  k=i  2  x  x  xk   j=o  1J   i=i  k=i  2  x   l0   lk   J=0  1J 

Letting  f  =  (-)0.p. , 

N   H  x    ,  k-1    . 

K  =  E   E  (-)0.^2I.n  -d   -2  E  d.  .[  , 
i=l  k=l  2  lL  l0   lk   j=0  ljJ 

and  noting  that 

H   k-1       H 
E  2  2  x. .  =   E  (2H-2k)x   , 
k=l   j=0  1J    k=l        lk 

the  total  cost  of  carrying  inventory  and  producing  the  N  products  over 
the  horizon  H  can  be  written  as 

N   H  N      H 

E   E  (2H-2k+lH\x.,  +  K  +   E  c   E  x  ,  (2.1) 

-ill  1  lk  1      ik 

i=l  k=l  i=l    k=l 

where  K  is  a  constant. 

The  system  controller's  objective  is  to  minimize  (2.1)  subject 
to  the  constraints  that  (a)  only  M  facilities  are  available  for  produc- 
tion during  each  period,  (b)  demands  must  be  met  without  allowing  back- 
orders,  and  (c)  at  most  one  product  can  be  scheduled  on  each  facility 
each  period  where  a  facility,  when  assigned  to  product  i,  is  assumed  to 
produce  p.  units  of  product  i.   A  mathematical  formulation  of  this 
problem  is  as  follows: 


15 


N  H  N      H 

Minimize   E  E  (2H-2k+l)i|f .  x.,  +  K  +  2  c   Ex 

i=l  k=l        x  lk     i=i  i  k=1  Ik 

subject  to 


[P2-1]  E  x   <  M  k  =  1,2,.. .,H 

i=l   lk 

k 

v      ^  l  =  1,2,. . . ,N 
E  x   S  w  '  '    ' 

.  1   ij     ik  k  =  1,2,. .. ,H 


x.,  :>  0,  integer   *  ~  1'2'"-  'N 
ik  6    k  =  1,2,... ,H 


It  is  easily  shown  that  problem  P2.1  has  a  feasible  solution  if  and  only 

N 

if   E  w   <  Mk   for  all  k  =  1,2, ,H. 

i=l   lk 

Note  in  problem  P2.1  that  all  costs  are  positive  and  that  any 

H 
feasible  solution  (X)  to  the  problem  must  have  Ex   s  w  '  for  all 

,  1   ik    iH 
k=l 

i  =  1,2,...,N.   Hence,  it  follows  that  any  optimal  solution  will  have 

each  product  i  produced  exactly  w   times  during  the  horizon.   This 

iH 

N     H 

implies  that  the  term,   E  c   E  x    is  a  constant  in  any  optimal  solu- 

i=  1   k= 1 

tion  to  problem  P2.1.   Knowing  that  constants  do  not  affect  the  optimi- 
zation, problem  P2.1  can  be  reformulated  as  follows: 

N   H 
Minimize   E   E  (2H-2k+l)\|r  x 
i=l  k=l 

subject  to 

N 
[P2.2]  E  x   <  M        k  =  1,2,...  ,H 

i=l 

k 

v      -.  i  =  1,2,.. .  ,N 

E  x  .  s  w  '  '    ' 

lj     ik      k  =  1,2,. . . ,H 

^  n  .    ,  i  =  1,2, . . . ,N 

x.  .  s  0,  integer        '  '    ' 
ik  6     k=l,2,...,H. 
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This  problem  can  be  viewed  as  a  flow  problem  which  can  be  solved 
using  any  of  the  algorithms  for  finding  minimal-cost  flows  in  single 
commodity  networks.   For  specific  details,  see  Ford  and  Fulkerson  [5], 
An  example  network  for  a  M-ma chine,  two-product,  three-period  problem 
is  shown  in  Figure  2.2.   Note  in  this  network  that  the  number  of  machines 
does  not  affect  the  number  of  nodes  and  arcs.   Rather  for  an  N-product , 
H-period  problem,  there  are  H(N+l)+2  nodes  and  H(2N+1)  arcs.   In  this 
network  the  variable  x   represents  the  flow  from  node  k  to  node  ik  along 
the  arc  (k,ik).   Efficient  algorithms  are  available  for  solving  this 
type  of  problem  (e.g.,  Fulkerson  [6]).   However,  because  of  its  special 
structure,  a  more  direct  solution  procedure  can  be  derived. 

A  Solution  Procedure 

The  problem  formulation  calls  for  each  product  i  to  be  assigned 
at  least  w   times  during  the  horizon  H,  and  since  all  costs  are  posi- 
tive, each  product  will  be  assigned  exactly  w   times  in  any  optimal 

in 

solution.   For  all  periods  the  increase  in  cost  for  assigning  product 
i  to  a  facility  in  period  k   as  opposed  to  a  facility  in  period  k  where 
k  <  k  <  H  is  2(k-k  ) ty . .   Hence,  from  the  standpoint  of  product  i,  it 
is  economically  more  advantageous  to  assign  product  i  in  period  k  rather 
than  period  k'. 

Suppose  that  the  products  are  numbered  so  that 

If  there  is  an  optimal  solution  to  problem  P2.2,  then  there  is  an 
optimal  solution  (X)  to  problem  P2.2  such  that 
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Lower 
Bounds 


Upper 
Bounds 


Figure  2.2   Network  Example  of  Problem  P2.1. 
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if  p  >  q  and 


if   Xpk<  ^pH—p.k-^-   E   Xpj  Where    Z  Xij^°'    (2'2> 

j=k+l   PJ       i=H+l   J 


then  x   =0. 
qk 

To  see  that  this  result  is  true,  assume  that  an  optimal  solution  (X*) 

to  problem  P2.2  has 

*  H    * 

x<(w-w)-£x. 

pk     PH    p,k-l       +1   pj 

and  x   >  0.   Let  k'  be  the  latest  period  such  that  k'  <  k  and  x   /  >  0. 
4K  pk 

Such  a  period  must  exist  if  the  solution  is  feasible.   By  making  a 
pairwise  swap,  a  feasible  solution  (X)  can  be  constructed  where 

for  (i=p,j=k)  and  (i=q,j=k') 

s    ,   .  .         for  (i=p,j-k')  and  (i=q,j=k) 

otherwise. 

This  solution  has  objective  function  value 

z  =   z     +   2(k-k')ilf      -   2(k-k')ilr 

q  TP 

* 
<  z      since      \|i      >  ilf    . 
P         q 

H 
This  process  can  be  continued  until  either  x   =  (w   -w     )  -   Ex 

pk    pH  p,k-l      +1  pj 

or  x   =  0.   In  either  case  an  optimal  solution  with  respect  to  p  and 
q  has  been  constructed.   Similarly,  by  starting  with  k  =  H  and  then 
letting  k  =  H-l ,  k  =  H-2,... ,k  =  2,  an  optimal  solution  can  be  con- 
structed such  that  condition  (2.2)  is  satisfied  for  all  products  and 
periods. 
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It  should  also  be  noted  that  once  the  products  are  numbered 
there  is  at  most  one  feasible  solution  which  satisfies  condition  (2.2). 
Hence,  if  a  solution  can  be  constructed  which  satisfies  condition  (2.2) 
for  all  products  and  periods,  the  solution  is  optimal. 

An  algorithm  based  on  the  above  result  can  now  be  given  for 
solving  problem  P2.2.   In  the  algorithm  each  planning  period  is  con- 
sidered once  where  period  H  is  the  first  period  assigned.   During  each 
period,  the  priority  in  which  assignments  of  products  to  facilities  are 
made  is  based  on  the  costs  iji..   The  larger  the  cost,  the  greater  the 
product's  priority.   Considering  each  product  in  the  order  of  its  pri- 
ority, the  number  of  facility  assignments  made  on  product  i  in  period 
k  is  based  on  (a)  the  availability  of  facilities  after  other  products 
with  a  higher  cost  have  been  assigned  and  (b)  the  maximum  desirable  num- 
ber of  assignments  of  product  i  in  period  k.   Formally,  the  algorithm  is 
as  follows  where  M  denotes  the  availability  of  facilities  at  intermedi- 
ate stages  of  the  procedure  and  p  and  t  indicate  respectively  the  cur- 
rent product  and  period  under  consideration. 

Step  0.   Order  the  product  indices  so  that 

ill      S:   ill  >    ...    2:   ill      ^   i|f    . 

rN         rN-l  *2         fl 

Define  w._  =  0  for  all   i  =  1,2 N,  and  set  t  =  H. 

iO 

Step  1.   Set  M'  =  M  and  p  =  N. 

Step  2.   Determine: 

H 
x   =  min  {m',  (w   -  w     )  -   E   x  .} 
Pt  pH     P,t-1     ,-_t+l   PJ 

M'  =  M'  -  x  . 
pt 
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Step  3.   Set  p  =  p-1.   If  p  S  1,  go  to  Step  2.   Otherwise, 

go  to  Step  4. 

Step  4.   Set  t  =  t-1.   If  t  H,  go  to  Step  1.   Otherwise, 

go  to  Step  5. 

H 
Step  5.   Terminate.   If   £  x.,  =  w   for  all  i  =  1,2,...,N, 

tel   lk    lH 
then  an  optimal  solution  (X)  to  problem  P2  has  been 

determined.   Otherwise,  the  problem  does  not  have 

a  feasible  solution. 

\ 
The  algorithm  is  a  one  pass  procedure  which  terminates  after 

exactly  HN  iterations  since  each  product  is  considered  once  during  each 
period.   At  termination,  the  constructed  sequence  is  deemed  to  be 
feasible  or  infeasible  depending  on  whether  or  not  w.   facility  assign- 
ments have  been  made  for  each  product.   Obviously,  an  infeasible  solu- 
tion can  be  made  feasible  only  if  the  number  of  facilities  is  increased 
or  if  the  no  stockout  restriction  is  relaxed.   Optimality  of  a  con- 
structed, feasible  solution  follows  from  the  result  given  in  (2.2). 

An  Example 

As  an  example  of  the  solution  procedure,  consider  a  two-product, 

three-period  problem  where 

M  =10  [2   3  12" 

W  = 
4  =  (10,20)  3   6   8 


The  product  indices  are  already  ordered  so  that   *   >  iji  .   Thus  the 
solution  procedure  goes  as  follows: 
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a.  Assign  product  2  twice  to  period  3  since  min  {10,2}  is  2. 

b.  Assign  product  1  eight  times  to  period  3  since  min  {8,9}  is  8. 

c.  Assign  product  2  three  times  to  period  2  since  min  {10,3}  is  3. 

d.  Assign  product  1  twice  to  period  2  since  min  {7,2}  is  2. 

e.  Assign  product  2  three  times  to  period  1  since  min  {10,3}  is  3. 

f.  Assign  product  1  twice  to  period  1  since  min  {7,2}  is  2. 

2  2 


;.   Stop.   An  optimal  solution  is   X  = 


function  value  z  =  760. 


3  3  2 


with  objective 


A  Warehouse  Space  Limitation 

In  many  production  environments,  there  exists  a  space  limitation 

on  the  maximum  amount  of  material  which  can  be  carried  in  inventory  at 

any  given  time.   In  the  problem  considered  in  this  chapter,  a  warehousing 

restriction  of  this  type  can  easily  be  incorporated  into  the  model  for 

the  special  case  where  one  machine-period  of  production  of  product  i 

takes  up  the  same  amount  of  space  as  one  machine-period  of  production 

of  any  other  product.   For  this  special  case,  let  S   denote  the  maximum 

k 

number  of  machine-periods  of  production  which  can  be  carried  in  inven- 
tory during  period  k  where  S   is  assumed  sufficient  to  store  the  initial 
product  inventories.   For  the  inventories  to  be  within  limits,  then, 
it  must  be  the  case  that  the  inventory  on  hand  at  the  end  of  each  period 

k  is  less  than  or  equal  to  S,  and  S,  1.   Equivalently ,  this  implies  that 

k  k+l 


N        k  N        k 

(2.3) 

M       :•«  .        '■ '  -'       i 

J=l      "  i=l   3- 


E  S  x.  .<min   {s,  ,S,    ,  }+     E       E   (d.  ,/p . )  -    E   (I .  _/p . ) 

.    ,  .    .    ij  k'    k+lJ       .  .  ij      i  lO      i 

i=l  j=l      °  i=l   .1=1        J  i=l 

for    all  k   =    1,2,. . . ,H 


22 

Obviously,  (2.3)  can  never  be  satisfied  if  the  right-hand  side 

is  strictly  negative.   Hence,  we  are  only  interested  in  pursuing  the 

problem  for  the  case  where  the  right-hand  side  is  nonnegative  for  all 

k  .=  1 ,2, .  .  .  ,H.   Also,  x.,  must  be  integer  valued  for  all  i  =  1,2 N 

lk 

and  k  =  1,2,...  ,H.   Therefore,  we  can  define  U   as  the  largest  integer 
less  than  or  equal  to  the  right-hand  side  of  (2.3).   The  inventory  con- 
straint can  then  be  included  in  problem  P2.2  by  annexing  the  constraint 

N   k 

£   E  x.  .  <  U,  . 
......  ij    k 

i=l  J=l 

The  revised  problem  is  still  a  network  problem  which  can  be 
solved  using  any  of  the  algorithms  for  finding  minimal-cost  flows  in 
single  commodity  networks.   In  fact,  this  warehousing  limitation  can 
easily  be  incorporated  into  the  previously  given  solution  procedure  by 
setting  X  =  0  in  step  0  and  in  step  1  redefining 

N   r 


M'  =  min  {M,     min    (U   -   S   Z   x  .)}. 
r=tf...  ,H   r    i=l  JUt 


Extensions  of  the  Model 

There  are  a  number  of  other  realistic  aspects  which  can  be 
easily  incorporated  into  the  model.   With  the  exception  of  (a)  below, 
the  previous  algorithm  is  not  satisfactory  for  determining  an  optimal 
solution  when  any  or  all  of  these  additional  enrichments  are  included 
in  the  model.   However,  an  optimal  solution  can  easily  be  determined 
for  these  cases  by  a  minimum-cost  flow  algorithm.   The  enrichments 
include: 
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(a)  The  number  of  facilities  available  in  each  period  may  vary. 
This  assumption  is  easily  included  in  the  previously  given  solution 
procedure  by  substituting  M   for  M  in  step  1  where  M  denotes  the 
number  of  facilities  available  in  period  t. 

(b)  In  addition  to  the  cost  of  producing  each  product  for  one 
machine-period,  there  can  be  a  cost  (convex)  incurred  on  the  number  of 
facilities  in  use  at  any  given  time.   This  can  be  introduced  into  the 
network  interpretation  of  the  problem  by  using  the  standard  trick  of 
replacing  each  arc  (s,k)  by  M  different  arcs  {(s,k)  •  I   =  1,2,...,M} 
each  having  an  upper  bound  of  one  on  the  flow  through  that  arc.   Any 
flow  on  arc  (s,k)  .  would  incur  a  cost  A.,  the  difference  in  cost  for 
using  I   facilities  during  period  k  rather  than  l-l   facilities.   Any 
cost  for  not  using  any  facilities  is  simply  a  constant  which  does  not 
affect  the  optimization.   The  convex  cost  assumption  implies  that 

AMSAM-1^  '•■  ^A2^A1 

which  insures  that  any  flow  in  period  k  will  be  placed  on  arc  (s,k). 
before  arc  (s,k)  .   unless  equality  holds  in  which  case  it  does  not 
matter. 

(c)  There  can  be  a  technological  restriction  on  the  number  of 
facilities  (m.)  which  can  be  used  to  pi'oduce  product  i  simultaneously. 
This  assumption  can  be  introduced  into  the  network  interpretation  of 
the  problem  by  placing  upper  bounds  of  m.  on  all  arcs  (k,ik). 

(d)  There  can  be  physical  limitations  (L.)  on  the  maximum  number 
of  units  of  product  i  which  can  be  carried  in  inventory  at  any  given 
time  so  long  as 
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k 

^v-nf*   „[Vik  +  !iO-   E  dij3     i  =  l,2,...,N. 
k=0,l,...,H  j=0   J 


This  assumption  can  be  introduced  into  the  network  interpretation  of  the 
problem  by  placing  upper  bounds  u   on  the  arcs  (ik;i,k+l)  where 


r     _     k       -, 

u.,  =  I  (L.  -  I   +   Id.  .)/p. 
ik   L  i    i0      ij  HiJ 


i  =  1,2,... ,N 
k  =  1,2,.. . ,H-1 


J 
and  [a]  is  defined  as  the  largest  integer  less  than  or  equal  to  a. 

(e)  Period  dependent  production  costs  (c   )  can  be  included  in 

ik 

the  network  interpretation  of  the  problem  by  adding  c.   to  the  cost 
(2H-2k+l)'i.  already  on  the  arcs  (k,ik). 

(f)  In  the  formulation  of  the  problem  P2.2,  it  has  been  assumed 
that  M  facilities  have  been  allocated  to  the  production  process  for  use 
throughout  the  entire  horizon  H.   Another  realistic  consideration  is 

a  problem  where  there  is  competition  for  these  facilities  on  a  long 
term  basis  (i.e.,  the  value  of  M  becomes  a  variable  in  the  optimization). 
This  competition  might  be  denoted  by  a  nonnegative  charge  f(M)  where 
logically  f(M)  would  be  a  monotonically  nondecreasing  function  of  M. 
When  this  charge  is  included,  a  procedure  for  problem  optimization  is 
to  solve  problem  P2.2    for     the  maximum  possible  value  of  M.   From 
the  optimal  solution  to  this  problem,  determine  the  maximum  number  of 
facility  assignments  (r)  in  any  period.   Then  an  optimal  solution  to  the 
allocation  problem  occurs  for  that  value  of  M  (M  =  0,1, . . . , r)  which  min- 
imizes z(M)  +  f(M)  where  z(M)  is  the  corresponding  optimal  solution  to 
problem  P2.2.   This  requires  that  problem  P2.2  be  solved  up  to  r  times. 

(g)   The  cost  of  producing  product  i  for  one  machine-period  can  be 
dependent  on  the  facility  as  well  as  the  period  in  which  the  assignment 
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is  made.   For  this  case  the  problem  formulation  becomes 

N   H  M        N   H   M 

Minimize   E   E  (2H-2k+l)i .   Ex.  .,  +  E   E   £  c  .  x  . 
i=l  k=l        x   jUl  l£k  1=1  k=l  jfci  Uk  iXk 

subject  to 

N   M 
[P2.3]  E   Ex    <M  k=l,2,...,H 

i=l  Jtl   k 


E   E  x.  .  .  s  w. 


1,2,...  ,N 


j=l  A=i  iij    ik       k  =  1'2>---  >H 

i  =  1,2,. ..  ,N 

XiXk  =  °  or  1  l  =    1,2,... ,M 

k  =  1,2,... ,H 

where  c  .   is  the  cost  of  producing  product  i  on  facility  I   in  period 
k  and  x  .  denotes  the  number  of  times  product  i  is  assigned  to  facil- 
ity I   in  period  k.   An  example  network  interpretation  of  this  problem 

is  given  in  Figure  2.3.   In  this  example  the  variable  x  .   represents 

iXk 

flow  from  node  m .,  to  node  n.,  along  the  arc  (m.  ,n   ). 
-ok  i  k  Xk   ik 

(h)  Realistically,  it  may  be  technologically  impossible  to  produce 
product  i  on  facility  Z.  This  restriction  can  be  included  in  the  model 
by  omitting  arcs  (m«k»nk)  in  aH  periods  k  =  1,2,...,H. 

Conclusions 

A  multi-facility,  production-inventory  scheduling  problem  has 
been  considered  which  can  be  solved  as  a  minimum-cost  flow  problem. 
However,  because  of  the  nature  of  the  objective  function,  an  even 
simpler  solution  procedure  has  been  developed  which,  even  for  large 
problems,  can  easily  be  carried  out  by  hand.   In  addition,  the  model 
can  be  generalized  to  allow  for  the  inclusion  of  other  realistic 
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Figure  2.3   Network  Example  with  Facility 
and  Period  Dependent  Costs. 
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assumptions.   Efficient  minimum-cost  flow  algorithms  are  available  for 
solving  the  general  problem  when  any  or  all  of  these  enrichments  are 
included. 

In  this  model  it  has  been  assumed  that  inventories  change 
linearly  from  the  end  of  one  period  to  the  end  of  the  following  period 
as  a  result  of  constant  production  and  demand  rates.   However,  the 
results  (with  the  exception  of  the  warehouse  space  limitation)  are 
equally  valid  for  any  system  in  which  inventory  changes  from  the  end 
of  one  period  to  the  end  of  the  following  period  do  not  drop  below  the 
minimum  of  the  beginning  and  ending  period  inventories.   This  includes 
the  inventory  assumptions  of  a  batch  type  production  system  in  which 
inventories  are  updated  only  at  the  end  of  each  production  period. 
The  only  revision  in  the  model  would  be  in  the  assessment  of  inventory 
carrying  charges. 


CHAPTER  3 


A  NETWORK  APPROACH  TO  THE  PRODUCT I ON- INVENTORY 
SCHEDULING  PROBLEM  WITH  BACKORDERING 


Introduction 

In  Chapter  2  the  multiple  facility,  multiple  product  production- 
inventory  scheduling  problem  without  backordering  was  considered.   In 
this  chapter  we  address  the  more  general  problem  where  backordering  to 
meet  demands  is  allowed.   Specifically,  the  cost  of  producing  each 
product  is  dependent  on  the  facility  on  which  the  product  is  produced 
as  well  as  the  period  during  which  the  assignment  is  made.   Production 
completed  during  period  k  is  added  to  each  product's  inventory  at  the 

end  of  the  period,  and  outstanding  backorders  (B.      units)  and  non- 

i ,  k— 1 

negative  demands  (d   units)  for  product  i  are  satisfied  from  this 
IK 

inventory.   Period  and  product  dependent  inventory  carrying  charges 
and  backordering  costs  respectively  are  assessed  on  material  remaining 
in  inventory  and  outstanding  demands  at  the  end  of  each  period.   During 
the  horizon  the  system  controller  must  schedule  the  products  to  the 
facilities  in  such  a  manner  as  to  meet  all  product  demands  on  a  first- 
come,  first-served  basis,  and  at  the  same  time  build  up  each  product's 

inventory  to  a  desired  nonnegative  ending  amount  (!.   units).   The  prob- 

iH 

lem  objective  is  to  determine  a  production  schedule  which  minimizes  the 
sum  of  production,  backordering  and  inventory  carrying  charges  over  the 

horizon. 
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In  succeeding  sections  we  formulate  this  multi-product,  multi- 
facility  production  scheduling  problem  as  a  linear,  mixed  integer  pro- 
gram.  Once  formulated,  this  problem  can  be  given  a  network  flow  inter- 
pretation which  looks  very  similar  to  the  network  interpretation  given 
by  Zangwill  [18]  for  his  one-facility  problem.   However,  in  contrast  to 
Zangwill's  problem,  this  mixed  integer  programming  formulation  can  be 
solved  as  a  minimum-cost  flow  problem  only  for  a  very  special  case. 
We  show,  however,  that  this  problem  can  be  reformulated  so  that  the 
general  case  can  also  be  solved  using  any  of  the  algorithms  for  finding 
minimal  cost  flows  in  single  commodity  networks.   This  reformulation, 
in  addition  to  being  a  model  for  a  very  general  multi-facility,  multi- 
product  problem  can  be  enriched  to  include  numerous  realistic  aspects 
of  the  production  system  without  severely  complicating  the. solution 
procedure. 

A  Mixed  Integer  Programming  Formulation 

Initially,  we  will  assume  that  there  is  a  nonnegative  number 

of  backorders  on  the  books  for  each  product  (i.e. ,  B   ^  0)  and  that 

initial  product  inventories  for  each  product  i  are  zero  (i.e.,  i".  =0). 

This  assumption  causes  no  loss  in  generality  since  the  first-come, 

first-served  rule  implies  that  the  first  i"   units  of  demand  for  each 

lO 

product  i  will  be  satisfied  from  the  initial  product  inventories  irre- 
spective of  any  decisions  made  in  the  problem  optimization.   Hence, 
we  can  assume  that  any  initial  inventories  have  been  used  to  adjust 
the  amount  of  initial  backorders  and  demands  appropriately  so  that 
initial  inventories  can  be  assumed  to  be  zero. 
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The  formulation  of  this  production  scheduling  problem  is  further 

facilitated  by  redefining  demands  in  terms  of  cumulative  machine-periods 

of  production.   This  is  done  by  defining  w   as  the  number  of  machine- 

ik 

periods  of  production  of  product  i  required  to  satisfy  all  demands  in 
periods  one  through  k.   Hence, 


k 

d „ 

ij     i0J  'i  k  =  1,2,... ,H-1 


ik     .  ,  ij     i0J  '  i 


j  =  l 


and 


{Ed   +  B,_  +  I.u}/P.      i  =  1,2,. ...N. 

,  ,  ik    lO    iHJ   i 


1H    k=l 


Define  x  as  the  number  of  times  product  i  is  assigned  to 

1JK 

facility  j  in  period  k,  and  denote  by  c. .,  the  cost  of  each  of  these 
assignments  where  c.    Sr  0.   Define  I    as  the  number  of  machine-periods 

IJ  K  IK 

of  production  of  product  i  in  inventory  at  the  end  of  period  k,  and 
define  B   as  the  amount  (in  machine-periods  of  production)  of  unsatis- 
fied demand  for  product  i  at  the  end  of  period  k  (i.e.,  the  total  amount 

of  product  i  on  backorder  at  the  end  of  period  k) .   Now  let  0        be  the 

ik 

cost  incurred  per  machine-period  of  production  of  product  i  in  inven- 
tory at  the  end  of  period  k,  and  let  b.,  be  the  cost  incurred  per 

ik 

machine-period  of  production  of  unsatisfied  demand  for  product  i  at  the 

end  of  period  k  where  0.,  >  0  and  b.,  >  0.   The  system  controller's 
ik  ik 

problem  is  to  minimize  the  sum  of  production,  inventory  and  backorder- 
ing  costs  over  the  planning  horizon  subject  to  the  constraints  that 
(a)  all  demands  must  eventually  be  satisfied,  (b)  only  one  product  can 
be  assigned  to  a  facility  each  period  where  a  facility,  when  assigned 
to  produce  product  i,  is  assumed  to  produce  exactly  p.  units  of  the 
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product,  and  (c)  the  inventory  level  for  each  product  i  at  the  end  of 

the  horizon  must  be  built  up  to  a  minimum  ending  amount  of  I   units 

iH 

A  mathematical  formulation  of  this  problem  is  as  follows: 

N   M   H  N  H-l         N  H-l 

Minimize   E   £   Z  c.  .,  x.   +  £   Ef,I  ,  +  E   E  b   B 

1=1     j  =  l    k=l     ljk     ljk        i„l     k=l     lR     lk        i=l    fal     ik     ik 


subject   to 


M       k 


S      2  x.  .     -  I      +  b.,   =  w  ,  i  =  J-.a.-'-iN        (3  1} 

j  =  l    4=1    1J  lk  ik  ik  k=    1,2,...,H-1      *•    *    ; 

M        H 

[P3-1]  ^        2  x.         ;>  w.„  i   =    1,2,...  ,N 

j  =  l    fel    ljk  lH 

N 

£  x..,    <  1  J   =   I.".-.- ,H 

i=1    ijk  k  =    1,2,.  .  .  ,H 

I       ,    B        >0  i    =    1.2.;..,N 

ik        ik  k  =    1,2,. . . ,H-1 

i  =   1,2,... ,N 

xiik   ^  °>    integer  j    =    1,2,. . . ,M 

k   =    1,2,... ,H 


Problem  P3.1  is  a  linear,  mixed  integer  program  which  can  be 

given  a  network  flow  interpretation.   To  construct  the  network,  define 

a  node  m   corresponding  to  each  facility  j  and  period  k,  and  define 

a  node  nik  corresponding  to  each  product  i  and  period  k.   From  each  node 

m   ,  construct  a  forward  arc  (m.,  ,n.,  )  to  each  node  n   for  all 
JK  jk   ik  ik 

i  =  1,2,...,N.   The  cost  per  unit  of  flow  on  each  of  these  arcs  is  c 

ijk 

Next  connect  each  node  n.,  to  node  n.  ,  ,  by  two  parallel  forward  arcs 

ik  i,k+l 

(nik>ni)k+1)  and  one  reverse  arc  (n    i'11^'   The  flow  through  one  of 

the  forward  arcs  must  be  exactly  w   .   Hence,  upper  and  lower  bounds  of 

ik 

Wik  are  Placed  on  the  flow  through  that  arc.   Each  unit  of  flow  on  the 
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other  forward  arc  incurs  an  inventory  charge  (0   ) ,  and  each  unit  of 

ik 

flow  on  the  reverse  arc  incurs  a  backordering  charge  (b   ).   Finallv 

ik  ' ' 

define  a  source  node  s  and  a  sink  node  t  where  a  forward  arc  (s  m  ) 

'  ik 

from  node  s  to  node  m   is  constructed  for  all  j  =  1,2,..., M  and 

JK 

k  =  1,2,... ,H  and  a  forward  arc  (n.T  t)  from  node  n   to  node  t  is  con- 

iH  iH 

structed  for  all  i  =  1,2,...,N.   Since  only  one  product  can  be  assigned 
to  a  facility  each  period,  a  capacity  of  one  is  placed  on  the  flow 

through  each  of  the  arcs  (s,m.,),  and  since  w.   machine-periods  of 

JK  iH 

product  i  must  be  produced  during  the  horizon,  a  lower  bound  of  w   is 

iH 

placed  on  the  flow  through  the  arc  (n..t).   Then  x    denotes  the  flow 

iH  ijk 

from  node  m   to  node  n   along  the  arc  (in   ,n.,),  I.,  and  w   corre- 
JK  ik  jk   ik    ik      ik 

spond  to  the  amount  of  flow  from  node  n   to  node  n      along  each  of 

ik  i,k+l 

the  forward  arcs  (n   ,n     ),  and  B.,  denotes  the  flow  from  node  n 

IK    1  ,  K+  J.  1 K  i  k_j_  1 

to  node  n   along  the  reverse  arc  (n.     ,n.,  ).   An  example  network  for 
•LK  i,k+l   ik 

a  two-product,  one-facility,  three-period  problem  is  given  in  Figure  3.1. 

In  this  problem  if  the  w   '  s  are  all  integer,  then  problem  P3.1 

can  be  solved  using  any  of  the  algorithms  for  finding  minimal-cost  flows 

in  single  commodity  networks.   For  specific  details  see  Ford  and 

Fulkerson  [5].   However,  because  of  the  way  in  which  the  w   's  have 

ik 

been  defined  (i.e.,  in  cumulative  machine-periods  of  production),  it  is 
unlikely  that  these  will  all  be  integers.   In  this  case,  using  a  minimal- 
cost  flow  algorithm  will  not  necessarily  provide  an  optimal  solution  to 
the  problem  which  has  all  x    integer.   In  the  remainder  of  this  chap- 
ter, we  will  show  how  problem  P3.1  can  be  reformulated  so  that  it  can  be 
solved  using  any  of  the  algorithms  for  finding  minimal -cost  flows  in 
single  commodity  networks  when  the  w   's  are  not  all  integer. 


33 


Production 
Costs 

Backordering 


Upper  Gnd 
Lower  Bounds 


nventory 
Costs 


Upper 
Bounds 


Figure  3.1   Network  Example  of  Problem  P3.1. 
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Reformulation  as  an  All  Integer  Program 

An  optimal  solution  to  problem  P3.1  has  the  property  that 
inventory  carrying  charges  and  backordering  costs  will  not  occur  simul- 
taneously for  any  product  i  and  period  k  (i.e. ,1   B   =0).   This 

ik  ik 

result  follows  from  the  fact  that  I  ,  and  B   occur  in  only  one  con- 

lk      ik  ' 

straint,  and  the  cost  coefficients  associated  with  each  variable  are 

strictly  positive. 

Consider  an  optimal  solution  to  problem  P3.1  having  B   >  0 

ik 

for  some  product  i  and  period  k.   Then,  from  the  above  result  I   =0 

ik 

and  from  constraint  (3.1) 

M   k 

E   E  x    +  B   =  w.  (3.2) 

j=l  Z=l   ljX    ik    lk 

Now,  breaking  up  each  term  into  its  integer  part  and  its  fractional 
part  where  [r]  denotes  the  largest  integer  less  than  or  equal  to  r  and 
f(r)  denotes  the  quantity  r  -  [r]  (i.e.,  for  nonnegative  r  the  fractional 
part  of  r  is  denoted  by  f(r)  where  0<  f(r)  <  1),  equation  (3.2)  becomes 

M   k 

E   E  x.    +  [B..  ]  -  [w.,  ]  =  f  (w*   )  -  f  (B.,  )  .  (3.3) 

._1  t   .  ljX     ik      ik       ik       ik  v    ' 

The  left-hand  side  of  equation  (3.3)  is  made  up  of  integer  terms 

only,  which  implies  that  f(w   )  -f  (B   )  must  be  an  integer.   This, 

i  k      i  k 

coupled  with  the  fact  that  0  <  f(w   )  <  1  and  0  <  f (B.  )  <  1,  implies 

ik  ik 

further  that  f(w   )  -f(B.,)  =  0.   Hence,  an  optimal  solution  to  problem 

ik      ik  * 

P3.1  having  B  ,  >  0  has 
ik 

f (B.,  )  =  f(w..  )  .  (3.4) 

ik       ik 
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In  similar  fashion  consider  next  an  optimal  solution  to  problem 

P3.1  having  I   >  0  for  some  product  i  and  period  k.   Then  B   =0  and 
ik  ik 

from  constraint  (3.1) 
M   k 

E   E  xiU  "  :ik  =  *ik'  (3-5) 

j=l  JLl   XJX/         lk    lk 


Breaking  up  each  term  into  its  integer  part  and  its  fractional  part, 
equation  (3.5)  becomes 

M   k 

.  ,  o-.    ij  *     lk      ik       ik       ik 

As  in  the  backordering  case,  the  left-hand  side  of  equation 

(3.6)  is  made  up  of  integer  terms  only  which  implies  that  f(w   )  +f(I   ) 

ik      ik 

must  be  an  integer.   This,  coupled  with  the  fact  that  (a)  f(w   )  equals 

ik 

zero  when  w   is  an  integer  and  is  strictly  positive  when  w    is  not 
ik  ik 

an  integer,  (b)  0<  f(w.,)  <  1,  and  (c)  0<  f (T   )  <  1,  implies  further 
ik  ik 


th 


at  f (w   )  +  f (I   )  equals  one  if  w    is  not  an  integer  and  equals  zero 
ik      ik  ik 


otherwise.   Hence,  an  optimal  solution  to  problem  P3.1  having  I   >  0 
has 

!1  -  f  (w   )    ,  if  w  t-    integer 
ik  ik 

(3.7) 
0  ,  otherwise. 

Finally,  by  defining  l"      as  a  0-1  variable  which  takes  on  the 
ik 

value  one  if  I   >  0  and  w   ^  integer  and  zero  otherwise,  relationships 
(3.4)  and  (3.7)  become 
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f  (B.,  )  =  f  (w   )  {l  -  I*  } 
lk       ik        ikJ 

and  (3.8) 

fdik)  =  fi-f(wik)}irk. 

Sufficient  relationships  have  now  been  developed  so  that  a  new 

formulation  of  the  production  scheduling  problem  can  be  given.   To 

formulate  this  problem,  let  I.,  and  b'   correspond  to  [I      1  and  TB  1 

lk      ik  ik      L  ik  ' 

respectively,  and  for  notational  convenience  define  a   =  f(w   )  and 

ik  ik 

Wik  =    [Wik]'      NoW   rePlacing    Iik   and   Bik    in  problem  P3.1   by   their   integer 
and  fractional   parts    and   then   substituting   the   expressions    in    (3.8)    for 
f^Bik^    and   f^I-k^'    WS  obtain   the   following   formulation: 

N        M        H  N     H-lf 

Min     Z       E       E  c        x        +    S       E10.  ,l'     +  { (1-a  ,  )0      -a     b      ll"} 
i=l  j=l  k=l   1J       Jk     i=i  k=l  i       ik     ik  ik     ikJ 


N     H-l 

+  E  E  b..  (b'  +  a.,) 
...,.,  ik  ik  ikJ 
i=l   k=l 


subject   to 


M        k 

S        E  X...-I'     -I*    +B'    =  w.  i   =   ^2,...,N 

1=1    &zl    1J  x  lk        ik        lk  k  =   1,2,..., H-l 


(3.9) 


M        H 
[p3-2]  E        Ex lik^w  i=l,2,...,N  (3.10) 

j  =  l   k=l    1Jk  lH 

N 

Ex..,    <     1  J  «   1.2,... ,11 

i=1    iJk  k  =  1,2,... ,H 

fl      ,    if    I        >  0   and  w.,   x  integer  i  =    1 ,2, .  .  .  ,N 

Iik  =  <  (3.11) 

lp      ,    otherwise  k  =    1,2,..., H-l 


Xijk'    Xik'    Bik  *  °>    integer  J    =    1,2,. ...M 


i  =  1 ,  2  , .  .  .  ,  N 
j  =  1,2,... ,M 
k  s    1,2,... ,H 
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Problem  P3.2  is  an  all  integer  program.   It  follows  from  our 
previously  derived  results  that  solving  problem  P3.2  is  equivalent  to 
solving  problem  P3.1.   We  now  proceed  to  show  that  problem  P3.2  can  be 
converted  into  a  minimal-cost  flow  problem. 

Formulation  as  a  Minimal -Cost  Flow  Problem 

Note  in  the  objective  function  of  problem  P3.2  that  the  term 

N  H-l 

T.       £  a   b    is  simply  a  constant,  and  since  constants  do  not  affect 
1=1  k=l 

the  optimization,  this  term  can  be  omitted.   Note  also  that  any  feasible 

solution  to  problem  P3.  2  must  have  all  x  .   integer.   Hence  w   in  con- 

ijk     °  ill 

straint  (3.10)  can  be  replaced  by  the  smallest  integer  greater  than  or 

equal  to  w    (i.e.,  (w   ))  without  any  effect  on  the  optimization, 
lri  In 

Finally,  we  will  show  that  constraint  (3.11)  can  be  replaced  by 


„  i  =  1  2     M 

I.,  =  0  or  e  ,  ,,...,«  (3.12) 

ik         ik        k  =  1,2,. . . ,H-1 


where  e.,  is  a  constant  defined  to  be  one  if  w   -  integer  and  is  zero 
ik  ik 

otherwise. 

To  see  that  (3.11)  can  be  replaced  by  (3.12),  let  P3.2'  denote 
problem  P3.2  with  (3.12)  used  in  place  of  (3.11).   From  equation  (3.5), 

I.,  >  0  if  and  only  if 

ik  ' 

M   k 

£   E  x.  .  ,  >  w.,  s  w.,  ,  (3.13) 

j  =  l  jtl  ljX    lk    lk 

and  from  constraint  (3.2)  and  the  integrality  requirements  of 

x.  .,  ,  I.,  =  0  if  and  only  if 
ljk   ik  J 
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jfl  |f1XiJ*SwikSwik"  (3'14) 

Together,  (3.13)  and  (3.14)  imply  that  constraint  (3.11)  is  equivalent 

to 


r 


M 


lk   i 


1   »   if   £   £  X.  .  .  -w.  >0  and  r,  ^  integer 
j  =  l  ^1  1J*   lk         lk 

i  =  1,2,. ..  ,N 
k  =  1,2,... ,H-1. 


,   otherwise 

(3.15) 

Note  next  that  problem  P3.2'  is  a  relaxed  version  of  problem  P3.2. 
Hence,  an  optimal  solution  to  problem  P3.2'  is  also  an  optimal  solution 
to  problem  P3.2  if  it  is  feasible  to  problem  P3.2.   From  (3.15)  this 
implies  that  an  optimal  solution  to  problem  P3.2'  is  also  an  optimal 
solution  to  problem  P3.2  if  in  this  solution 

M   k 

(a)  l"      ■=   1   implies   £   £  x.  .  .-w.,  >0  and  w   ^integer 

ik  j=1  jtl  iji   lk        ik 

and 

M   k 

(b)  i''   =  0   implies   £   £  x.  .  -w.,  <0  and/or  w   =  integer. 

ik  j  =  1  jtl  iji   ik  ik 

Case  (a) : 

Obviously,  w   cannot  be  integer  valued  since  this  would  imply 

that  I    would  be  restricted  to  the  value  zero.   Hence,  to  prove  this 

„  M   k 

case  we  assume  that  I   =  1,  w.,  ~  integer,  and   £   £  x    -w   <0  for 

lk       lk  j=l  fcl  ^Z        ik 

some  product  i  and  period  k  in  an  optimal  solution  to  problem  P3.2'. 

Then  in  order  for  constraint  (3.9)  to  be  satisfied,  it  must  be  the  case 

that  Bik  >  1.   This,  however,  implies  that  there  exists  a  feasible 
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solution  to  problem  P3.2   having  l"      and  b'   each  reduced  in  value  bv 

ik      ik  J 

one  yielding  an  objective  function  value  improvement  given  by 

Az  =  -  (i-a.,)0.1  +  a.b,  -  b., 
ik  *ik    ik  ik    ik 

<  0   ,   since  (1  -  a.,  )  ,  b.,  ,  0   >  0. 
ik    ik    ik 

Hence  a  contradiction  since  we  could  not  have  started  with  an  optimal 
solution. 

Case  (b) : 

Obviously,  (b)  holds  if  w   =  integer.   Hence  to  prove  this  case 

M   k 
we  assume  that   i"  =  0,  w  4  integer,  and   £   Z  x   „  -  w   >  0  for  some 
lk      lk  j=l  jtl  ljX   ik 

product  i  and  period  k  in  an  optimal  solution  to  problem  P3.2'.   Then  in 

order  for  constraint  (3.9)  to  be  satisfied,  it  must  be  the  case  that 

Iik  s  1.   This,  however,  implies  that  there  exists  a  feasible  solution 

to  problem  P3.2   having  I    increased  in  value  by  one  and  i'   decreased 

ik  ik 

in  value  by  one  yielding  an  objective  function  value  improvement  given 

by 

Az  =  -0    +  (1  -a   )0    -  a.,  b.. 
ik        ik   ik    ik  ik 

<  0   ,   since  a.,  ,  b   ,0   >  0. 
ik'   ik   *ik 

Hence  a  contradiction  since  we  could  not  have  started  with  an  optimal 
solution. 

Therefore,  cases  (a)  and  (b)  do  hold  at  optimality  in  problem 
P3.2  which  implies  that  solving  problem  P3.2'  is  equivalent  to  solving 
problem  P3.2.   Making  the  specified  changes  in  the  objective  function 
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and  constraint  (3.10)  in  problem  P3.2  and  replacing  constraint  (3.11) 
by  (3.12),  problem  P3.2  becomes  as  follows: 

N       M       H  N     H-lf 

Min       E       2       E  c.  .,  x.  .,    +      E       E  10.,  i! 


i=l   j=l   k=l    iJk   ijk        i=l   fc=lL    ik   ik 


f^-ik^i^iAk^lJ  +    A    ^hiA 


N      H-l 

-i        R 

ik 


i=l   k=l 


subject   to 

M       k 


i   =    1 ,  2  ,..'.,  N 


.Z,     //ijU  '   Iik"Iik+Bik"Wik  k  =  1,2 H-l 

j=l    *=1 

M        H 

[P3.3]                         ^      Z^*  <w.H)  1=  1,2,...  ,N 

I  j=  1,2,. ..,M 

.    ..    iJk  k  =  1,2,. . . ,H 
i^l 

z«     <e  i  =  1,2 N 

ik           ik                                       '  k  =  1,2,. . . ,H-1 


x..,,    I.,,    I.,,    B,     SO,    integer  j    =    1,2,.  .  .  ,M 

ljk        ik        ik        ik 


i  =  1,2,. .. ,N 
j  =  1,2,... ,M 
k  =  1,2,. . . ,H 


The  multiple  facility,  multiple  product  production-inventory 
scheduling  problem  with  backordering  is  now  in  the  format  of  problems 
which  can  be  solved  using  any  of  the  algorithms  for  finding  minimal- 
cost  flows  in  single  commodity  networks.   Efficient  algorithms  are 
available  for  solving  this  type  of  problem  (e.g. ,  Fulkerson  [  6  ]). 
The  network  structure  corresponding  to  problem  P3.3  is  similar  to  that 
of  problem  P3.1.   The  primary  difference  is  in  the  construction  of 

a  second  arc  (n.,  ,n.  ,  .,  )  from  each  node  n  ,  to  node  n  .  „  .   Flow 
ik   i,k+l  ik  i,k+l 

through  each  of  these  arcs  is  limited  by  the  capacity  e   ,  and  a  cost, 
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{ (1-a  ,  )0.  -  a  b   ] ,  is  incurred  for  each  unit  of  flow  placed  on  the 

1  J'i    J.  K     JL  K  l.K 

arc  (n   ,n.     ) .   The  actual  amount  of  flow  placed  on  this  arc  is 
denoted  by  the  variable  I*.   An  N-product,  M-facility,  H-period  problem 
of  this  type  will  have  H(MfN)+2  nodes  and  H(M+MN+4N) -3N  arcs.   An  example 
network  is  shown  in  Figure  3.2  for  a  one-facility,  two-product,  three- 
period  problem  where  for  notational  convenience  s    is  defined  to  be 

ik 

f(1-aik)0ik-aikbik5- 

In  addition  to  being  a  model  for  a  very  general  multi-facility, 

multi-product  production-scheduling  problem,  there  are  a  number  of  other 
realistic  enrichments  to  the  problem  which  can  be  included  in  the  model 
without  severely  complicating  the  solution  procedure.   These  are  dis- 
cussed below: 

(a)  There  can  be  a  technological  restriction  on  the  number  of 
facilities  (q.,  )  which  can  be  used  to  produce  product  i  simultaneously 
in  period  k.   This  generalization  can  be  introduced  into  problem  P3.3 
by  adding  the  constraint 

M 

v      -  i  =  1 , 2 , .  .  .  ,  N 

ljk    ik  k  =  1,2,. . .  ,H. 

(b)  In  addition  to  the  cost  of  producing  each  product  for  one 
machine-period,  there  can  be  a  cost  incurred  on  the  number  of  facil- 
ities in  use  at  any  given  time.   Let  r  denote  the  cost  per  facility 

k 

in  use  during  period  k  where  r  >  0.   Then  this  enrichment  is  included 

k 

in  the  problem  formulation  by  adding  the  following  term  to  the  objec- 
tive function, 

H    N    M 
Er,  E    £  x.  M  . 
k=l  k  i=l  j=l  ljk 
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Production  Upper  and  Upper 

Costs  Lower  Bounds     Bounds 

nyentory 

^  posts 

Adjustments 

Costs 
e, 


Upper 
Bounds 


Figure  3.2   Network  Example  of  Problem  P3. 3. 
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However,  this  cost  function  does  not  have  to  be  linear.   Rather  the 
costs  can  be  convex.   The  insertion  of  a  convex  cost  term  is  a  stan- 
dard trick  in  networks  and  is  discussed  in  detail  in  Ford  and  Fulkerson 
[  5]. 

(c)  The  number  of  facilities  available  in  each  period  may  vary. 
This  generalization  is  incorporated  in  the  model  by  omitting  those 
nodes  m    (and  the  associated  arcs)  in  period  k  corresponding  to  any 
facility  j  which  is  unavailable. 

(d)  Upper  bounds  can  be  placed  on  the  values  of  i'   and  B   implv- 

lk      ik 

ing  physical  limitations  on  the  amount  of  product  i  which  can  be  carried 
in  inventory  or  on  backorder  during  period  k+1. 

(e)  Realistically,  it  may  be  technologically  impossible  to  produce 

product  i  on   facility  j.   This  restriction  can  be  included  in  the  model 

by  omitting  arcs  (m  ,  ,n. ,  )  in  all  periods   k  =  1 , 2, . . .  ,H. 
jk   ik 

Conclusions 

A  multiple  facility,  multiple  product  production-inventory 
scheduling  problem  with  backordering  has  been  considered.   In  a 
straightforward  fashion  this  problem  has  been  formulated  as  a  linear, 
mixed  integer  program  which  can  be  given  a  network  flow  interpretation. 
However,  this  formulation  can  be  solved  as  a  minimal-cost  flow  problem 
only  for  a  very  special  case.   In  this  paper  it  has  been  shown  that  the 
general  problem  can  be  reformulated  so  that  it  can  be  solved  using  any 
of  the  algorithms  for  finding  minimal-cost  flows  in  single  commodity 
networks.   This  reformulation,  in  addition  to  being  a  model  for  a  very 
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general  multi-facility,  multi-product  problem,  can  be  enriched  to  include 
numerous  other  realistic  aspects  of  the  production  system  without 
severely  complicating  the  solution  procedure. 


CHAPTER  4 


A  MULTIPLE  FACILITY,  MULTIPLE  PRODUCT  PRODUCTION 
SCHEDULING  PROBLEM  WITH  OVERTIME 


Introduction 

In  the  real  environment,  many  production  systems  face  fluctuating 
demand  patterns  which  often  tax  the  productive  capability  of  the  system 
to  meet  these  demands.   As  a  consequence,  many  schedulers,  or  system  con- 
trollers, use  overtime  to  increase  the  productive  capacity  of  the  system. 
When  this  option  is  available,  the  cost  of  the  overtime  must  be  taken 
into  account,  and  production  costs  become  a  function  of  the  amount  of 
overtime  used  as  well  as  the  straight  time  production  amounts.   In  this 
chapter  the  problem  of  scheduling  lots  to  facilities  is  considered  for 
the  case  where  an  overtime  option  is  available  to  the  system  controller. 

Specifically,  an  industrial  process  made  up  of  M  facilities  in 
parallel  is  considered  where  there  are  N  different  products  to  be  pro- 
duced over  a  finite  planning  horizon.   During  any  period  k,  M  facilities 
are  available  for  production  where  M  <  M.   All  of  these  facilities  are 
considered  to  be  identical  in  that  p.  units  per  period  of  product  i  can 
be  produced  by  any  available  facility.   During  any  given  period  at  most 
one  product  can  be  scheduled  on  each  facility.   Corresponding  to  each 

of  these  assignments  is  a  cost,  denoted  by  c. .,  ,  which  includes  the  cost 

ijk 

of  using  facility  j  in  period  k  to  produce  product  i  as  well  as  the 
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associated  costs  of  producing  p.  units  of  product  i.   At  the  end  of 
each  production  period,  the  system  controller  has  the  additional  option 
of  scheduling  an  overtime  shift  for  producing  product  i  on  facility  j 
when  product  i  is  scheduled  on  straight  time  on  facility  j.   The  incre- 
mental cost  per  period  for  this  overtime  production  is  denoted  by  g. 
where  g^^  >  0.   (Alternatively,  this  problem  can  be  thought  of  in  a  batch 
type  production  process  where  the  scheduler  must  decide  whether  to  pro- 
duce zero,  one,  or  two  batches  of  product  each  time  the  product  is 
assigned  to  a  facility  for  one  period.   It  can  also  be  viewed  in  the 
context  of  a  continuous  production  process  where  the  scheduler  has  the 
capability  of  doubling  the  production  rate  of  each  of  the  facilities.) 
Production  completed  on  straight  time  and  overtime  during  period  k  is 
added  to  each  product's  inventory  at  the  end  of  the  period,  and  nonneg- 
ative  demands  (d   units)  are  satisfied  from  this  inventory.   During 
the  horizon  the  system  controller  must  schedule  the  products  to  the 
facilities  in  such  a  manner  as  to  meet  all  product  demands,  without 
suffering  backorders  or  stockouts,  and  at  the  same  time  build  up  each 
product's  ending  inventory  level  to  a  minimum  nonnegative  amount  denoted 
by  IiH>   The  problem  objective  is  to  determine  such  a  production  sched- 
ule which  minimizes  the  sum  of  straight  time  and  overtime  production 
costs  over  the  horizon. 

In  succeeding  sections,  the  overtime  production  scheduling 
problem  is  formulated  as  a  linear,  integer  program.   This  formulation 
can  be  given  a  network  flow  interpretation,  but  in  general  the  problem 
cannot  be  solved  using  any  of  the  minimal-cost  flow  algorithms.   How- 
ever, it  is  shown  that  this  problem  can  be  solved  by  solving  a  relaxed 
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version  of  the  problem.   This  relaxed  problem  has  the  property  that  it 
can  be  solved  using  any  of  the  very  efficient  algorithms  for  finding 
minimal-cost  flows  in  single  commodity  networks.   Finally,  a  single  pass 
algorithm  is  developed  for  solving  a  special  case  of  the  overtime 
problem. 

The  Overtime  Production  Scheduling  Problem 

In  this  problem  product  inventories  are  updated  only  at  the  end 
of  each  production  period,  and  product  demands  must  be  satisfied  from 
these  end-of-period  inventories  without  allowing  backordering  or  stock- 
outs.   Hence,  requiring  that  the  net  ending  period  inventories  be  non- 
negative  for  all  products  and  periods  is  a  necessary  and  sufficient  con- 
dition for  all  demands  to  be  met  without  backorders  or  stockouts. 

To  facilitate  the  formulation  of  this  problem,  define  w   as 

ik 

the  minimum  number  of  machine-periods  (integer)  during  which  product  i 
must  be  produced  in  the  first  k  periods  in  order  to  have  a  nonnegative 
inventory  level  at  the  end  of  period  k.   Then, 

k 
w   =  max  [0,  <(  £  d..-7   )/p)]        i=l,2,...,N 
lk  j=1   iJ     l0    i  j        k  =  1,2,. . . ,H-1 

where  <a>  denotes  the  smallest  integer  greater  than  or  equal  to  a. 
In  addition  to  being  nonnegative,  product  i's  inventory  level  at  the 
end  of  the  horizon  must  be  greater  than  or  equal  to  I   .   Hence,  define 


ill' 

H 

r'ij  "  'iO-V 


w.h=  max  [0,  <(i   +   E  d   -  l.Q)/p.)}  i  =  1,2,...,* 


j  =  l 

and  let  W  be  an  X  by  H  matrix  with  components  Tw   1 

L  ikJ 
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Define  y.  .,  as  the  number  of  times  facility  j  is  scheduled  to 
ijk  J 

produce  product  i  on  straight  time  during  period  k,  and  denote  by  c 

ijk 

the  cost  of  each  of  these  assignments  where  c    >  0.   Similarly, 

ijk 

define  x.    as  the  number  of  times  facility  j  is  scheduled  to  produce 

IJK 

product  i  on  overtime  during  period  k,  and  denote  g.  as  the  incremental 
cost  (product  dependent  only)  of  producing  product  i  for  any  overtime 
machine-period  where  g.  >  0.   Then  the  system  controller's  problem  is 
to  minimize  the  sum  of  straight  time  and  overtime  production  costs  sub- 
ject to  the  constraints  that  (a)  only  M  facilities  are  available  for 
production  in  period  k,  (b)  demands  must  be  satisfied  without  allowing 
backorders  or  stockouts,  (c)  product  i  can  be  scheduled  on  overtime  on 
facility  j  during  period  k  only  if  it  has  already  been  scheduled  on 
straight  time  on  facility  j  during  period  k,  and  (d)  only  one  product 
can  be  scheduled  on  a  facility  each  period  where  a  facility,  when 
assigned  to  produce  product  i,  is  assumed  to  produce  exactly  p.  units 
of  the  product  each  machine-period.   A  mathematical  formulation  of  this 

problem  is  as  follows: 

N   H   Mk 
Minimize   E   E   E   f  c .  .,  y .  .,  +  g .  x.   "} 
1=1  k=l  j  =  l    ljk   ljk     X  ljk' 


subject  to 


jB=1  j=l   1J*    1J* 


ik 


[P4.1] 


N 

2  y 
i=l 


ijk 


<  1 


ijk'   ijk 


y     >  x 

*ijk    ijk 


5  0,  integer 


i   =    1,2,. 

.  ,N 

k   =    1,2,. 

•  ,H 

k   =    1,2,. 

•  ,H 

j    -    1,2,. 

••"k 

i   =    1,2,  . 

•  ,N 

k   =    1,2,. 

•  ,H 

j    -    1,2,. 

•'\ 

i    =    1,2,  . 

•  ,N 

k  =    1,2,. 

•  ,H 

j    =    1,2,.. 

■>\ 

(4.1) 


(4.2) 


(4.3) 


(4.4) 
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Problem  P4.1  is  a  linear,  integer  program  which  can  be  given 
a  network  flow  interpretation.   However,  the  problem  cannot  in  general 
be  solved  using  any  of  the  minimal-cost  flow  algorithms.   We  will  show, 
though,  that  a  relaxed  version  of  the  problem  can  be  formulated  and 
that  this  relaxed  problem,  when  solved,  will  provide  an  optimal  solu- 
tion to  problem  P4.1.   This  relaxed  problem  will  have  the  property  that 
it  can  be  solved  using  any  of  the  algorithms  for  finding  minimal-cost 
flows  in  single  commodity  networks.   The  formulation  of  this  problem  is 
based,  in  part,  on  the  following  result: 

Lemma  4 . 1 

A  necessary  and  sufficient  condition  for  problem  P4.1  to  have 

a  feasible  solution  is  that  there  exists  a  set  of  y    's  satisfying 

ijk 

constraints  (4.2)  and  (4.4)  and  having 
k   Mi 

L  jfi  y^i>-  <wik/2>  <4-5> 

for  all  i  =  1,2,. . . ,N  and  k  =  1,2, ... ,H. 

Proof  of  Lemma  4.1 

A.   Necessary  Condition: 

Assume  that  there  exists  a  feasible  solution  (Y*,X*)  to 
problem  P4.1  having 

k     I 

E   Z  yiii  <  <Wik/2>  (4-6> 

1=1   j=l  X3Z  lk 

for  at  least  one  product  i  and  period  k.   Then  (4.6)  and  constraint 
(4.3)  together  imply  that 
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k     I 
4_1  j_l    XJ x    iJ-6       lk  ik 

A  contradiction  exists  since  constraint  (4.1)  cannot  be  satisfied  unless 
(4.5)  holds. 

B.   Sufficient  Condition: 

Set  x.    =  y.    for  all  i  =  1,2,. . . ,N;  k  =  1,2,. . . ,H;  and 
j  =  1,2, ...,11.   Then,  by  construction,  constraints  (4.1),  (4.3)  and 
(4.4)  are  satisfied.   Hence,  (Y,X)  is  a  feasible  solution  to  problem 
P4.1. 

Q.E.D. 

Since  condition  (4.5)  must  be  satisfied  by  all  feasible  solu- 
tions to  problem  Pl.l,  it  can  be  included  in  the  problem  as  a  redundant 
constraint.   To  formulate  a  relaxed  version  of  problem  P4.1,  then,  we 

add  constraint  (4.5)  to  the  problem,  omit  constraint  (4.3),  and  make 

H   Mk 
the  variable  transformation   z.  =   £   £  x.   .   The  relaxed  problem 

1   k=l  j  =  l  ** 
can  be  formulated  as  follows: 

N       H        \  N 

Minimize        £        £        £c..y  +£erz 

1=1   k=l   j=l      ljk      lJk        i=l      i    i 

subject   to 

k        M£ 

zi +  ^  .^ *m* wik  i:l:l:::::l    (4-7> 

[p4-2]  Ji^-1  j  =  i;22;:::;^    (4-8) 

M,, 
v        v  '  i    =    1,2,. . .  ,N  ,,    „„ 


■i 
Jfcl    J  =  l    "iji   S    <Wik/2>  k=    1'2'"--.H 


y .  z.    >   0,    integer  k   =    1,2,. . .  ,H 

ljk         l 


i  =  1,2,.. . ,N 
k  =  1,2,. . . ,H 
J   =   1,2,...,^ 
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Since  (a)  problem  P4.2  is  a  relaxed  version  of  problem  P4.1,  and 
(b)  problem  P4.2  includes,  as  constraints,  the  necessary  and  sufficient 
conditions  specified  by  Lemma  4.1,  it  follows  that  problem  P4.1  has  a 
feasible  solution  if  and  only  if  problem  P4.2  has  a  feasible  solution. 
Let  (Y,Z)  be  an  optimal  solution  to  problem  P4.2.   Then  (Y°,X°)  is  an 
optimal  solution  to  problem  P4.1  where,  by  construction, 


yijk  =  y 


ijk 


ijk 


<°   -  < 


Ll 


if 


k   j 

Z   E  y^ 
1=1   r=l 


irX 


<  z 


i   =    1,2,. 

.  ,N 

k   =    1,2,. 

. ,H      (4.10) 

J    =    1,2,. 

-\ 

otherwise. 


To  show  that  (Y  ,X  )  is  a  feasible  solution  to  problem  P4.1, 
note  first  that  constraints  (4.2),  (4.3)  and  (4.4)  are  satisfied  by 
construction  since  (Y,Z)  is  a  feasible  solution  to  problem  P4.2.   To 
show  that  constraint  (4.1)  is  satisfied,  we  note  from  (4.10)  that 
Y  =  Y  and  that  either 


k   ;"X        k   mi 
fcl  j=l  13i       fcl  j=l  lji 


(4.11) 


(b) 


k   "A 

E   E  x° 
A=l  j=l  1J 


(4.12) 


for  all   i  =  1,2,.  .'.  ,N  and  k  =  1,2,.  .  .  ,H. 
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Case  (a) : 

After  substituting  Y   for  Y,  constraint  (4.9)  and  equation  (4.11) 
together  imply  that 


Hence,  constraint  (4.1)  is  satisfied  for  this  case. 

Case  (b) : 

Constraint  (4.7)  and  equation  (4.12),  together  with  the  fact 
that  Y  =  Y  implies  that 

k  Mx   r 

£   E  {y .  .  .  +  x .  . }  s  w 
Z=l   j=l  lji  lj£     lk 

Hence,  constraint  (4.1)  is  also  satisfied  for  this  case,  and  (Y°,X°)  is 
a  feasible  solution  to  problem  P4.1. 

In  minimization  problems,  it  is  well  known  that  the  optimal 
objective  function  value  for  a  relaxed  version  of  a  problem  represents 
a  lower  bound  on  the  value  of  the  objective  function  of  any  feasible 
solution  to  the  unrelaxed  problem.   Therefore,  since  (a)  (Y°,X°)  is 
a  feasible  solution  to  problem  P4.1,  and  (b)  the  objective  function 
value  of  the  solution  (Y  ,X  )  is  equal  to  the  objective  function  value 
of  the  optimal  solution  (Y,Z),  it  follows  that  (Y°,X°)  is  an  optimal 
solution  to  problem  P4.1. 

Problem  P4.2  can  be  viewed  as  a  flow  problem  which  can  be  solved 

using  any  of  the  algorithms  for  finding  minimal-cost  flows  in  single 

commodity  networks.   An  example  network  for  a  one-product,  three-period 

problem  with  two  facilities  available  in  each  period  is  shown  in 

Figure  4.1.   In  this  network  the  variable  y.  .,  represents  the  flow  from 

ijk 
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node  m .,  to  node  n.,  along  the  arc  (m   ,n   ),  and  the  variable  z   rep- 
jk  lk  jk   ik  i 

resents  the  flow  from  node  s  to  node  o    along  the  arc   (s,o   ). 
Efficient  algorithms  are  available  for  "solving  this  type  of  problem 
(e.g. ,  Fulkerson  [  6 ]). 


A  Special  Case  of  the  Overtime 
Scheduling  Problem 


For  the  special  case  where  the  production  costs  are  not  facility 

and  period  dependent  (i.e.,  c..,  =  c.  for  all  k-  1,2,. ...H  and 

ijk    i 

j  =  1,2,... ,M.  ) ,  a  single  pass  procedure  can  be  developed  for  solving 
the  overtime  production  scheduling  problem.  This  procedure,  even  for 
large  problems,  can  easily  be  carried  out  by  hand.   Before  describing 

the  procedure,  however,  note  first  that  facility  discrimination  is  no 

"k 

longer  necessary,  and  we  can  let  y   =   E  y    .   Problem  P4.1,  then, 

.    1      ijk 

reduces    to    the   following   formulation: 

N        H  N 

Minimize        E        £     c.y.,    +      E     g   z 
i=l   k=l  i=l 


subject   to 

k 

■7         _i_  Y     t: 

iX  ik  k  =   1,2 H 


z-    +      £  7.  «>-  w-,  ?    ~~    }'„""  '?  (4.13) 


*        X=l 

[P4.3]  N 

Z  yik  ~  \  k=l,2,...,H  (4.14) 

i=l 

£yn  *  <wu/a>  l:l:l:::::l     (4-i5) 

^   _       .     .  i   =    1,2,. . . ,N 

y.,  ,    z.    so,    integer  '     ' 

iki  '  6  k=l,2,...,H 
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The  procedure  to  be  given  for  solving  problem  P4.3  involves 
two  phases.   Phase  I  determines  a  set  of  y   '  s  (if  they  exist)  which 
satisfy  the  necessary  and  sufficient  conditions  for  a  feasible  solution 
to  problem  P4.3  to  exist  as  stated  in  Lemma  4.1.   Phase  II,  then,  builds 
onto  this  Phase  I  solution  in  order  to  determine  an  optimal  solution 
to  problem  P4. 3. 

Problem  P4.3  has  a  feasible  solution  if  and  only  if  there  exists 
a  set  of  yik's  which  satisfy  constraints  (4.14)  and  (4.15)  for  all 

i  =  1,2,... ,N  and  k  =  1,2,. . . ,H.   The  set  of  y   '  s  is  determined  by 

lk  ' 

omitting  constraint  (4.13)  and  solving  the  relaxed  problem.   In  terms 
of  the  network  given  in  Figure  4.1,  this  is  equivalent  to  omitting  those 
nodes  o±k  for  all  i  =  1,2,. ..,N  and  k  =  1,2,. . . ,H  and  all  arcs  leading 
into  and  out  of  these  nodes.   To  solve  this  relaxed  problem,  a  backward 
procedure  is  used  which  is  essentially  equivalent  to  a  flow  algorithm. 
The  procedure  is  as  follows  where  R  denotes  the  availability  of  facil- 
ities at  intermediate  stages  of  the  procedure  and  p  and  t  indicate 
respectively  the  current  product  and  period  under  consideration. 


Phase  I 


Step  0.   Define  (wiQ/2>  =  0  for  all  i  =  1,2,... ,N  and  set  t  =  H. 
Step  1.   Set  R  =  M  and  p  =  N. 
Step  2.   Determine: 

ypt=  min  [R,  «wpH/2>-<wp(t_l/2»-   I  y   } 


j  =  t+l 


H 

where     E  y  .  =  0 


j  =  H+l  PJ 


R  =  R  -  y 

pt 
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Step  3.   Set  p  =  p-1.   If  p  >  1,  go  to  Step  2.   Otherwise, 

go  to  Step  4. 

Step  4.   Set  t  =  t-1.   If  t  >  1,  go  to  Step  1.   Otherwise, 

go  to  Step  5. 

H 
Step  5.   Terminate.   If   I  y   =  <w  /2>  for  all  i  =  1,2,..., N, 

k=l 

then  an  optimal  solution  has  been  determined  for  the 
relaxed  problem.   Otherwise,  problem  P4.3  does  not 
have  a  feasible  solution. 

If,  in  Phase  I,  an  optimal  solution  to  the  relaxed  problem  is 
determined,  then  we  proceed  to  Phase  II,  since  from  Lemma  4.1  a  feasible 
solution  to  problem  PI. 3  is  known  to  exist.   Otherwise,  we  stop,  since 
problem  P4.3  does  not  have  a  feasible  solution.   Before  proceeding, 
however,  note  that  any  feasible  solution  to  the  relaxed  problem  having 
each  product  i  assigned  exactly  <w   "2>  times  during  the  horizon  is  an 
optimal  solution  to  the  relaxed  problem.   Hence,  the  relaxed  problem 
may  have  many  alternate  optimal  solutions.   Using  the  backward  proce- 
dure approach  to  determine  one  of  these  optimal  solutions  has  the  fol- 
lowing important  advantage,  however,  in  determining  an  optimal  solution. 
The  Phase  I  procedure  at  termination  for  each  period  k  has  the  property 
that  either  (a)  all  v     facilities  are  assigned,  or  (b)  the  only  product 
assignments  which  can  be  reassigned  to  period  k  (without  violating 
feasibility  or  assigning  a  product  more  than  (w   ,2/  times)  are  those 
which  are  already  assigned  to  a  period  r  where  r  >  k. 
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Example  1 

As  an  example  of  the  Phase  I  procedure  for  solving  the  relaxed 
version  of  problem  P4.3,  consider  a  two-product,  three-period  problem 
where 


c  =  (6,3) 
M  =  (8,8,10) 


W  zr 


4    6   24 
6   12   14 


<W/2>   = 


2  3   12 

3  6    7 


An  optimal  solution  is  determined  as  follows: 

a.  Assign  product  two  once  to  period  3  since  min  {l,10}  =  1. 

b.  Assign  product  one  nine  times  to  period  three  since 
min  {9,9}  =  9. 

c.  Assign  product  two  three  times  to  period  two  since 
min  {3,8}  =  3. 

d.  Assign  product  one  once  to  period  two  since  min  [l,5]  =  1. 

e.  Assign  product  two  three  times  to  period  one  since 
min  [3,8]  =  3. 

f.  Assign  product  one  twice  to  period  one  since  min  {2,5}  =  2. 

"2  1  9~~ 


Stop.   An  optimal  solution  is   Y  = 


3  3  1 


Assuming  Phase  I  yields  a  feasible  solution  to  the  relaxed 
problem,  we  now  proceed  to  Phase  II.   In  this  phase  an  optimal  solution 
to  problem  P4.3  is  determined  by  building  on  to  the  Phase  I  solution. 
This  is  equivalent  to  assuming  that  the  part  of  an  optimal  solution 
which  has  been  determined  in  Phase  I  is  fixed  (i.e.,  let  y'  denote  the 
solution  obtained  in  Phase  I)  and  now  all  that  remains  is  the  determin- 
ation of  the  remaining  portion  of  the  solution  (i.e.,  denote  the  remain- 
ing portion  of  the  solution  by  (y',Z)).   The  Phase  II  problem  to  be 
solved  is  formulated  as  follows: 
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N   H  N 

Minimize   S   E   {c.y!   +  c.y'  }  +   Eg 
i=l  k=l    1  lk    x  lk    i=l  i 


subject  to 


z.+   £  y"   ;>  w   -   £y'  i=l,2,...,N 

1    4=1  U    ik    A=l  iX  k  =  1,2,. ...H 


[P4.4] 


■  Vik-  \  "  .E/ik  k=  1-2,. ...H 

1=1  1=1 

y*  ,  z.  2:  0,  integer  1  =  i.S,....  ,N 

lk   x  k  =  1,2,. ..,H 

A  network  interpretation  of  problem  P4.4  is  obtained  from 

Figure  4.1  by  omitting  all  arcs  (n.,  ,n.  ,  )    for  all  i  =  1,2 N  and 

ik   i ,k+l  '  '    '" 

k  =  1,2,. . . ,H-1,  and  all  arcs  (n.^.t)  for  all  i  =  1,2 N.   Also 

In  ' 

N 

redefine  \  =   \   -      £  yj[      for   all    k  =    1,2,...,H.      Now  note   from   this 

i=l 

network  interpretation  that  for  all  products  having  c.  >  g.  that  it  is 
at  least  as  advantageous  to  assign  overtime  as  it  is  to  assign  straight 
time.   Since  this  can  be  done  without  causing  facility  conflicts,  all 
demands  for  these  products  are  satisfied  by  overtime  assignments.   For 
those  products  having  c±  <   gi ,  it  is  more  advantageous  to  assign  straight 
time  than  it  is  overtime.   However,  there  are  a  limited  number  of  avail- 
able facilities  left.   Hence,  it  is  more  advantageous  to  assign  a  product 
i  having  (c^-gj  <  ^C,S.)    to  a  vacant  facility  than  it  is  a  product  j. 
To  determine  an  optimal  solution  to  problem  P4.4,  we  again  use  a  back- 
ward procedure  because  of  the  useful  property  which  it  allowed  in  the 
Phase  I  procedure.   In  the  procedure,  sets  S  and  S  are  defined  where  S 
contains  those  products  for  which  an  optimal  form  is  known  (i.e.,  ini- 
tially those  products  having  c   >  g.),  and  S  contains  those  products 
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for  which  an  optimal  solution  form  is  not  known.   The  procedure  then 

considers  one  product  in  S  at  a  timed,  e.,  the  first  product  considered 

is  that  product  in  S  having  the  minimum  value  of  c  -g  ) .   Once  the  form 

i   i 

of  an  optimal  solution  is  determined  for  this  product,  then  it  is  trans- 
ferred from  set  S  to  S.  This  process  continues  until  S  is  an  empty  set. 
In  the  procedure  we  build  on  to  the  y  's  directly  rather  than  using  the 
Y   and  Y   notation. 

Phase  II 

Step  0.   (a)  Order  the  product  indices  so  that 

Cl"gl  S  Vg2  *  •••  "  W 

(b)  Define  the  sets 

S  =  [if  c  -g  <  0}  and  S  =  {ijc.-g.  >  0}. 

H 

(c)  Set  z  =  w   -  E  y   for  all  i  =  1,2... . ,N. 

i    ui   k=1  ik 

Step  1.   If  S  ■£   0,  go  to  Step  2.   Otherwise,  stop..   An  optimal 

solution  is  (Y,Z). 

Step  2.   Set  p  =  max  i  and  t  =  H. 
¥i£S 

Step  3.   (a)  Determine 

N  t-1 

a  =  min  [z  ,  M  -   Ey.z  +   Zy.-w     ] 

where  yiQ  and  w    are  defined  to  be  zero  for  all 
i  =  1 , 2  , .  .  .  ,  N. 
(b)  Set 

y  4.  =  y  4.  +   a 
pt  •'pt 

and 
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Step  4.   If  t  =  1  or  z.  =  0,  go  to  Step  5.   Otherwise,  set 

t  =  t-1  and  return  to  Step  3. 
Step  5.   Set  S  =  {s-p}  and  S  =  [s+p],  and  go  to  Step  1. 

Since  overtime  periods  are  assigned  without  facility  conflict, 
it  is  clear  that  the  only  way  there  could  exist  a  better  optimal  solu- 
tion to  problem  P4.3  than  that  obtained  by  Phases  I  and  II  is  if  there 
exists  an  optimal  solution  (Y*,Z)  having 

H   *     H  H        H 

k=l      k=l  k=l  J    k=l 

for  at  least  two  products  i  and  j  where  c.-g.  <  c.-g..   However  such 
a  solution  cannot  exist  since  (a)  in  Phase  II  product  i  was  assigned 
to  all  periods  before  product  j,  and  (b)  in  Phase  I  the  backward  proce- 
dure has  the  property  that  there  does  not  exist  a  feasible  interchange 
of  these  two  product  assignments. 

Example  2 

As  an  example  of  the  Phase  II  procedure,  reconsider  the  problem 
given  in  Example  1  where,  in  addition,  g  =  (5,9).   An  optimal  solution 
is  determined  as  follows: 

a.  The  products  are  already  indexed  so  that  c.  -g   >  c  -e  . 

1  bl    2   2 

Hence,  since  c  -g  =  1  and  c^-g  =  -6,  we  define 
S  =  [l]  and  S  =  {2}. 

b.  Set  y   =1+0=1  since  rain  {7,0,1}  =  0. 

c.  Set  y   =3+4=7  since  min  {7,4,4}  =  4. 

d.  Set  y   =3+3=6  since  min  {4,3,3}  =  3. 


e.   Stop.   An  optimal  solution  is  Y 


2  1  9~| 
_6  7  lj 


and  Z  =  (12,0). 
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Conclusions 

A  multiple  facility,  multiple  product  overtime  scheduling 
problem  has  been  considered.   A  linear,  integer  programming  formulation 
was  given  for  the  problem  which  in  this  case  can  be  solved  by  solving 
a  relaxed  version  of  the  problem.   This  relaxed  problem  was  shown  to 
have  the  property  that  it  can  be  solved  using  any  of  the  algorithms  for 
finding  minimal-cost  flows  in  single  commodity  networks.   However,  for 
an  important  special  case  a  single  pass  algorithm  was  developed  which 
is  more  efficient  than  the  standard  flow  algorithms. 


CHAPTER  5 


AN  EFFICIENT  INTEGER  PROGRAMMING  ALGORITHM 
FOR  A  MULTIPLE  BATCH  PRODUCTION  SCHEDULING  PROBLEM 


Introduction 


In  this  chapter  a  multiple  batch  production  process  is  considered 
where  the  horizon  length  is  a  variable  in  the  optimization.   During  any 
production  period,  each  of  the  M  facilities  has  the  capability  of  pro- 
ducing up  to  mi  batches  of  product.   A  learning  effect  is  experienced 
in  the  production  of  these  batches  and  is  exhibited  in  the  production 
costs.   The  cost  (c  ^)  of  the  I        batch  produced  on  a  facility  during 
any  period  is  no  less  than  the  cost  (c..  ,   )  of  the  iflSt  batch  produced 
for  all  I   =  1,2, . . . ,m  -1.   A  batch  of  product  i  contains  p.  units,  and 
only  whole  batches  are  produced.   Production  completed  during  period  k 
is  added  to  each  product's  inventory  at  the  end  of  the  period,  and  non- 
negative  demands  (d   units)  for  product  i  are  satisfied  from  this 
inventory. 

It  is  assumed  here  that  under  normal  operating  conditions  the 
scheduler,  or  system  controller,  for  this  process  has  a  production  plan 
for  scheduling  the  X  products  to  the  M  facilities.   This  production  plan 
may  be  a  schedule  based  on  the  well-known  infinite  horizon  economic  manu- 
facturing quantity  which  was  discussed  in  Chapter  1.   in  these  formula- 
tions of  the  problem  all  demands  must  be  met  without  allowing  backorder- 
ing,  and  the  problem  objective  is  to  determine  a  repetitive  production 
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cycle  for  allocating  the  production  time  of  a  facility  to  each  of  the 
N  products  so  that  costs  are  minimized.   In  most  treatments  the  addi- 
tional assumption  is  made  that  initial  product  inventory  levels  are 
"in  phase  with"  the  infinite  horizon  cycle  plan  meaning  that  the  product 
inventory  levels  are  all  sufficiently  large  for  an  infinite  horizon 
cycle  plan  to  begin  immediately  without  suffering  stockouts.   When  ini- 
tial product  inventories  are  considered  in  the  problem,  there  may  exist 
a  subhorizon  of  unknown  time  duration  (h)  during  which  the  facilities 
must  be  scheduled  "out  of  phase  with"  the  infinite  horizon  cycle  plan 
in  order  to  avoid  stockouts  or  backorders.   In  this  chapter  an  extremely 
efficient  procedure  is  developed  for  bringing  a  system  into  phase  with 
a  given  production  pattern  such  as  an  infinite  horizon  cycle  plan. 

Let  IiH  donote  the  minimal  inventory  level  in  units  for  each 
product  i  which,  when  taken  collectively  for  all  products,  represents 
an  acceptable  condition  or  "state  of  the  system"  for  the  system  con- 
troller's desired  production  pattern  to  begin.   Realistically,  there  is 
usually  more  than  one  such  state  which  can  be  used  since  the  entry  point 
into  a  repetitive  production  cycle  is  ordinarily  not  unique.   (However, 
the  number  of  possible  entry  states  is  finite  since  time  is  discrete.) 
In  this  chapter  we  develop  a  solution  procedure  for  a  single  entry 
state  and  leave  to  the  reader  the  straightforward  reapplication  of  the 
solution  procedure  for  other  entry  states.   The  scheduler's  problem 
then  is  to  determine  an  assignment  of  products  to  facilities  and  the 
corresponding  production  amounts  at  each  of  these  assignments  over  a 
variable  horizon  h  subject  to  the  constraints  that  (a)  all  demands  are 
satisfied  without  allowing  backordering,  and  (b)  the  nonnegative, 
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initial  product  inventory  level  I    (units  of  product  i)  for  all 
products  is  converted  to  a  desired  nonnegative,  ending  inventory  level 
IiH»   The  objective  is  to  minimize  the  sum  of  production  costs  and  a 
general  cost  function  f(h)  which  is  a  monotonically  nondecreasing  func- 
tion of  the  horizon  length  h. 

In  succeeding  sections,  we  identify  a  specially  structured 
integer  program  and  present  an  extremely  efficient  algorithm  for  its 
solution.   After  identifying  special  properties  of  this  general  problem, 
the  multiple  facility,  multiple  product  production  scheduling  problem 
is  formulated  as  an  integer  program  and  shown  to  have  the  same  special 
structure. 

An  Integer  Program 

Noteworthy  contributions  have  been  made  to  the  present  state-of- 
the-art  of  solving  discrete  optimization  problems  (i.e.,  for  a  survey 
of  many  of  these,  see  Geoffrion  and  Marsten  [7  1),    but  to  date  the 
research  has  not  provided  an  efficient  procedure  for  solving  the  general 
problem.   Rather,  researchers  have  been  most  successful  in  (a)  identify- 
ing problem  structures  which  insure  integer  termination  of  the  corre- 
sponding relaxed  (non-integer)  optimization  problem,  and  (b)  determin- 
ing solution  techniques  which  are  applicable  to  specific  classes  of 
discrete  problems. 

In  this  section  a  specially  structured  integer  programming 
formulation  is  identified  for  which  a  very  efficient  algorithm  will  be 
presented  for  determining  its  solution.   The  problem  formulation  is 
given  below  where  (a)  J  is  the  index  set  of  integers,  (b)  all  functions 
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are  real  valued,  (c)  gk  (•)  exists  for  all  k  =  1 , 2, . . . ,m,  and  (d)  foi 

any  vectors  e  =  (e  ,e  ,... ,e  )  and  x  =  (x. , x  . . . . ,x  )  having  e   :>  0 

a  £•  n  —     12      n  i 

and  x±   >   L  for  all  i  =  1 , 2 , .  .  .  , n ,  f  (x)  <  f  (x  +  e)  and 


g~1{hk(x)}<g-1{hk(x+   e)}. 


Minimize   z  =  f(x) 
subject  to 


[P5.1]  g  (x.  )  -  h  (x)  >  0  k=l,2,...,m 

k  Jk     k  "  \  ^  J 


(5.1) 


x.  >  L. ,  integer        i  =  1,2, . . . ,n      (5.2) 

A  result  which  provides  the  basis  for  an  algorithm  for  solving 
problem  P5.1  is  the  following: 

Lemma  5.1 

If  lower  bounds,  L.  ,  are  known  for  all  variables  x   i  =  1,2 n. 

1  i       '  '    '  ' 

in  problem  P5.1,  then  for  the  variable  x.   in  any  constraint  k  of  (5.1), 

Jk 

\    -   max{L    (g^1  {hk(L)}>}  ,   jk  €  j  (5.3) 


is  also  a  lower  bound  on  the  variable  x.   in  problem  P5.1.   (Note  that 

Jk 
(a)  implies  the  smallest  integer  greater  than  or  equal  to  a.) 


Proof  of  Lemma  5.1 

Assume  that  there  exists  an  x .   <  L .  in  a  feasible  solution  to 

Jk    Jk 
problem  P5.1.   By  assumption,  L.   is  a  known  lower  bound  on  the  value 

Jk 
of  x   in  any  feasible  solution  to  the  problem.   Hence  it  must  be  the 

Jk 
case  that  x  >  L.   which  implies  from  (5.3)  that 
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L*  =  (g"1  [h  (L)}>  .  (5.4) 

Jk     k     k  - 

But  (5.4)  is  just  constraint  k  of  (5.1)  rewritten  as  an  equality  with 

each  variable  x.  set  at  its  known  lower  bound.   Since  g   fh  (L) ]  < 
i  k   k  —  J 

—1  r  -,  * 

g,  [h  (L  +  e) }  for  any  e  >  0,  the  only  way  a  variable  x.   <  L.   can  be 
k   k  —   — •  —   —  in 

Jk    Jk 

feasible  to  problem  P5.1  is  if  some  x  =  L   is  decreased  in  value.   Hence, 

l     l 

a  contradiction  exists  since  L.  is  a  lower  bound  on  the  value  of  x   in 

1  i 

any  feasible  solution  to  problem  P5.1. 

Q.E.D. 

An  Algorithm  for  Solving  the  Integer  Program 

As  is,  the  problem  formulation  has  an  interesting  structure. 
The  objective  function  is  a  nondecreasing  function  for  all  increasing 
values  of  x. .   Each  of  the  constraints  in  (5.1)  has  one  variable  which 
can  be  expressed  as  a  nondecreasing  function  for  all  increasing  values 
of  x..   A  feasible  solution  to  the  problem  requires  that  all  variables 
be  integer  valued  and  bounded  from  below.   Using  these  properties  and 
assuming  that  a  feasible  solution  to  the  problem  is  known  to  exist,  an 
algorithm  for  solving  the  problem  can  be  described  as  follows: 

Step  0.   Set  x  =  L. 

Step  1.   If  any  constraint  in  (5.1)  is  not  satisfied,  go  to 

Step  2.   Otherwise,  stop.   An  optimal  solution  is  (x). 

Step  2.   Pick  any  unsatisfied  constraint  in  (5.1),  redefine 

\  =  <s;>k(*)}> , 

and  return  to  Step  1. 
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At  each  iteration  of  the  procedure,  the  algorithm  either 
terminates  or  the  value  of  a  variable  in  an  unsatisfied  constraint  is 
increased  by  at  least  one.   Therefore,  since  (a)  a  feasible  solution 
is  known  to  exist,  and  (b)  from  Lemma  5.1  each  variable  is  always  reset 
at  a  value  equivalent  to  a  known  lower  bound  on  that  variable,  the 
algorithm  must  terminate  in  a  finite  number  of  iterations.   Furthermore, 
since  the  value  of  the  objective  function  in  problem  P5.1  is  nondecreas- 
ing  for  an  increasing  value  of  any  x. ,  optimality  of  the  solution  at 
termination  also  follows  from  Lemma  5.1. 

An  Example 

Consider  the  following  example  which  is  in  the  format  of 
problem  P5. 1. 

Minimize   z  =  0.5x  +  2x 

subject  to 

2X1  "  °-5X2  S  10  (5-5) 

"  Xl  +    4X2  ^   8  (5-6) 

x,  -  O.lx  >.      1 
1        2 

xi  '  xr>   S  °>  integer 
The  problem  is  illustrated  in  Figure  5.1  where  an  optimal  solution  is 
obtained  by  the  iterative  scheme  in  the  following  manner. 

a.  Set  x  =  0. 

b.  Constraint  (5.5)  is  not  satisfied  so  redefine 
xx  =  <{10  +  0.5(0)}/2)  =  5. 

c.  Constraint  (5.4)  is  not  satisfied  so  redefine 
x2  =  <{8  +  1(5) }/4)  =  4. 
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2x,-0.lx2  =10 


-x,  +  4x2=8 


Figure  5.1   Geometric  Interpretation  of  the  Algorithmic  Exampl. 
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d.  Constraint  (5.5)  is  not  satisfied  so  redefine 
K.    =  <{10  +  0.5(4)}/2>  =  6. 

e.  Stop.   All  constraints  are  satisfied.   An  optimal  solution 
is  x  =  (6,4)  with  objective  function  value 

z  =  0.5(6)  +  2(4)  =  11. 

Other  Properties  of  the  Integer  Program 

In  addition  to  the  ease  of  solution,  problem  P5.1  has  a  number 
of  additional  properties  which  are  of  interest. 

Property  1 

Consider  any  two  problems  having  identical  constraint  sets  of 
the  form  required  in  problem  P5.1,  but  different  objective  functions. 
The  same  optimal  algorithmic  solution  solves  both  problems  irrespective 
of  the  objective  function  f(x)  so  long  as  each  satisfies  the  condition 
f  (x)  <  f  (x  +  e)  for  all  x  and  any  e  >.   0. 

Property  2 

Any  feasible  solution  (x)  to  a  problem  p' ,  where  problem  P 
includes  the  constraint  set  of  problem  P5.1,  lias  x  >  x*  where  x*  is 
the  optimal  algorithmic  solution  to  problem  P5.1. 

Property  3 

The  algorithmic  solution  to  problem  P5.1  is  the  unique  optimal 
solution  to  the  problem  if  f  (x)  <  f (x  +  e)  for  all  x  and  all  e  ^  0 
where  at  least  one  component  of  e  is  strictly  positive. 


70 

Property  4 

The  corresponding  continuous  portion  of  problem  P5.1  {omitting 
the  integer  requirements  on  x  in  constraint  (5.2)}  has  an  optimal  solu- 
tion (x  )  which  has  all  x.  integer  when  a  feasible  solution  exists  if 
gk  [(\(*)1    is  inteSer  valued  for  all  x  =  integer  and  k  =  1,2,.  .  .  ,m. 

Note  finally  that  in  the  development  t    the  algorithmic  proce- 
dure, termination  was  insured  since  a  feasii  .  solution  to  problem  P5.1 
is  known  to  exist.   If  a  feasible  solution  is  not  known  to  exist,  then 
an  additional  stopping  rule  must  be  added  to  the  procedu  e.   This  can 
be  done  in  a  straightforward  manner  when  necessary  (e.g.,  place  upper 
bounds  on  each  of  the  decision  variables  and  terminate  the  solution  pro- 
cedure when  any  one  of  these  bounds  is  exceeded). 

The  Production  Scheduling  Problem 


The  system  controller's  problem  is  to  determine  an  assignment 

of  N  products  to  M  facilities  and  the  corresponding  production  amounts 

at  each  of  these  assignments  over  a  variable  horizon  h.   This  is  to  be 

done  in  such  a  manner  as  to  [a]  satisfy  all  demands  (d   units  for 

ik 

product  i  in  period  k)  without  allowing  backordering,  [b]  convert  the 

nonnegative,  initial  inventory  level  (I.Q)  for  all  products  i  into 

a  desired  nonnegative,  minimum  ending  inventory  level  (T.  ),  and  [c] 

iH 

minimize  costs.   The  costs  include  production  charges  (i.e.,  c.  denotes 

l  jL 

the  incremental  cost  of  producing  an  t        batch  of  product  i  during  any 
machine-period)  and  a  general  cost  function  f(h)  which  is  a  monotonically 
nondecreasing  function  of  the  horizon  length  h.   A  batch  of  product  i 
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contains  exactly  p.  units  of  the  product  and  only  whole  batches  can  be 

produced.   Only  one  product  can  be  assigned  to  a  facility  each  period, 

and  up  to  m  batches  of  that  product  i  can  be  produced  during  any 
i 

machine-period. 

By  assumption,  the  initial  and  required  ending  period  inventory 
levels  are  nonnegative  which  implies  that  backordering,  if  it  occurs, 
must  occur  at  the  end  of  at  least  one  of  the  intermediate  periods 
(k  =  1 ,2, . . .  ,h-l) .   Define  w   as  the  minimum  number  (integer)  of 
batches  of  product  i  which  must  be  produced  during  the  first  k  periods 
in  order  to  prevent  backorders  from  occurring  at  the  end  of  period  k. 
Then, 

wik  = max  [o> «  eAj  "Tio)/px^        i:  i,1:: :::*-!. 

Likewise,  define 


w..  =  max  {0,  ((I   +   Ed..  -  I.n)/p.»     i  =  1,2,...,N. 

lh  ih    .  ,  ij     iO    l  J 

Now  by  defining  x.  .   as  the  number  of  facilities  scheduled  to 
produce  exactly  JL   batches  of  product  i  during  period  k,  a  mathematical 
formulation  of  this  problem  ca n  be  given  as  follows: 


N   h   m±      I 
Minimize  f(h)  +   E   E   E   E  c.  x.  , 
i=l  k=l  1=1   r=l 


subject  to 


j=l  £=1 

[P5.2]  N   m. 


1  —  1  9.  W 

E   E   Xx.    >  w       X  -  ,*;••••  »J      (5.7) 
iXj     lk    k  =  1 ,2  , .  .  .  ,h 


E   Ex.,,   <  M      k  =  1,2,.  .  .  ,h      (5.8) 
i=l  X=l 


h ,  x.  .   SO,  integer      i  -   1 , 2 , . . .  ,mi 


i  =  1 ,  2  , .  .  .  ,  N 
I  =  1,2,. . . ,mj 
k  =  1,2,. . . ,h 
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Problem  P5.2  is  a  nonlinear,  integer  program  since  the  index  h 
is  a  variable  in  the  optimization.   Rather  than  solve  problem  P5.2 
directly,  a  relaxed  version  of  the  problem  will  be  formulated  and 
solved.   This  relaxed  problem  can  be  put  into  the  format  of  problem 
P5.1  and  solved  very  efficiently  by  the  previously  given  algorithm. 
Once  this  relaxed  problem  has  been  solved,  a  simple  procedure  for  solv- 
ing problem  P5.2  will  be  given. 

The  relaxed  problem  is  formulated  by  (a)  omitting  those  con- 
straints in  (5.7)  corresponding  to  k  =  l,2,...,h-l,  (b)  relaxing  the 

k  constraints  given  in  (5.8)  to  the  single  constraint  given  in  (5.10) 

h 

below,  and  (c)  making  the  variable  transformation  y  .  =   Ex., 

1 1       ,  ,  i  £k 
k=l 

The  problem,  then,  is  as  follows: 

m. 

N         i  I 

Minimize      f  (h)    +      E        £  l     c        y.  „ 

■19-1  i      ir  J1Z 

i=l    JL=1  r=l 

subject    to 


[P5.3]  £      Xy..>w  i=l,2,...,N  (5.9) 

4=1        lX  lh 

- (i™  jx  i  "u 


h,    y±Jt  ,  0,    xntcger        ^  __   ^         ' 


Theorem   5. 1 


A  sufficient  condition  for  problem  P5.3  to  have  a  feasible 

solution  is  that 

N 


£   (d.,  /m.p.)  <  M  (5.11)  . 

.  .,    lk   l  l 
i=l 


for  all  k  =  1,2,. . . ,h. 
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Relation  (5.11)  is  a  facility  capacity  condition  which  is 
common  in  continuous  infinite  horizon  production-inventory  problems. 
For  example,  see  Hadley  and  Whitin  [8  ].   A  proof  for  this  multiple 
facility,  multiple  batch  discrete  problem  is  given  in  the  Appendix 
where  the  proof  given  is  along  the  lines  of  that  given  by  Elmaghraby 
and  Mallik  [2  ]  for  the  case  m.  =  1  and  M  =  1.   In  this  paper  it  is 
assumed  that  condition  (5.11)  is  satisfied. 

Note  that  problem  P5.3  is  in  the  format  of  problem  P5.1  for 

the  special  case  where  m.  equals  one  for  all  i  =  1,2,...,N.   However, 

in  general,  it  is  not  in  the  correct  format.   Define  the  variable  v  . 

±JL 

as  the  number  of  machine-periods  having  at  least  £   batches  of  product 
i  produced.   Then,  for  all  i  =  1,2,...,N, 

\l  ~   Vi,£fl     ,  I   =  1,2,..  .  ,m4-l 

i 

'ii 

Making  this  variable  transformation  in  problem  P5.3,  we  obtain  the 

following  equivalent  program. 

N   mi 
Minimize   f(h)  +   E   E  c  .v  . 
i-1  Jtl  U   ±l 

subject  to 

mi 

E  v  .  s  w  i  =  1,2,...  ,N 

j6=1  lZ         lh 

[P5.4]  v    _  v       s  0  i  =  1,2,. ..,N 

ii    l,i+l  I  =1,2,.. . ,m.-l 


h  :>  (1/M)  S  v.  , 
i=l 


.        -^  n   -  *  1=1, 2, ...,N 

h,  v.  .  2:  0,  integer  .     '  '    ' 

ii  i  =  1,2,. .. ,m.  . 
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!*   __.,  13*  I   =  1,2,.  .  .  ,m.  , 


Now  making  the  additional  transformation, 

£ 
r=l 

problem  P5.4  can  be  reformulated  equivalently  as  follows: 

N  l 

Minimize   f(h)  +   £  [c.    t.     -   £   (c   .   -  c  ,)t  ,1 
i=l   1,mi  X'mi    JL=1        1>'e+1     i*   iJtJ 


subj  ect  to 


t.     s  w  i  =  12 N 


[p5-5^  2t±1  -  t12  £  0  i  =  1,2 N 


t. .  - 1 


i  =  1,2,.. . ,N 


I  i,A-l  i  =    2,3,. 


N 
h  >  (1/M)  £  t. 
i=l  x 

h ,  t .  .  >  0 ,  integer 


1 ,  2 , .  .  .  ,  N 


iX     '     s  X  =  1,2,. 

X 

Note  that  for  the  case  m.  =  2,  constraints  (5.12)  should  be  omitted. 

The  constraint  set  of  problem  P5.5  is  now  in  the  format  of 
problem  P5.1  (i.e.,  each  constraint  has  one  variable  which  can  be 
written  as  a  nondecreasing  function  of  other  nondecreasing  variables 
and  a  feasible  solution  is  known  to  exist).   Hence,  an  optimal  solution 
to  problem  P5.5  can  be  obtained  by  the  algorithmic  procedure  since 

c.1>c.0^...>c     >-  0 
il    i2  i,m 

l 

implies  that  the  objective  function  is  a  nondecreasing  function  of  any 
increasing  variable. 
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Upon  solving  problem  P5.5  by  the  algorithmic  procedure,  the 

optimal  solution  (h,T)  can  be  transformed  back  into  an  optimal  solution 

(h,Y)  to  problem  P5.3.   This  solution  has  the  property  that  m  batches 

i 

of  each  product  i  will  be  produced  each  time  product  i  is  assigned  for 

one  machine-period  with  the  possible  exception  of  one  assignment. 

m  -1 

1 

Hence,      S     y..   equals   0   or   1    for   all    products    i  =    1,2,...  ,N. 
jfc=l 


Lemma   5.2 


The  optimal  algorithmic  solution  (h,T)  to  problem  P5.3  has 
m  -1 


E  y  ,  =  0  or  1  (5.13) 

A=l    X 

for  all  products  i  =  1,2,...,N. 

Proof  of  Lemma  5.2 

From  Property  1,  the  optimal  algorithmic  solution  to  problem 

P5.3  is  the  same  irrespective  of  whether  the  objective  function  value 

is  strictly  increasing  or  just  nondecreasing  in  each  variable.   Hence, 

in  problem  P5.3  it  can  be  assumed  that 

c...  >  c   >  .  ..>  c 
ll     i2         m. 

l 

for  all  products  i  =  1,2,...,N.   Assume  that  the  optimal  solution  to 
problem  P5.3  does  not  satisfy  (5.13).   Then  an  optimal  solution  (h*,Y*) 

must  have   E  y.  .  >   2  which  implies  that  either 

JL=1      lZ 

m.  -1 

1     * 

(a)  £  iy.  .   =  9  <  m. 

jfcl    Xi 

or 

m  -1 
1     * 

(b)  I  £y.      =  6  >.   m.  . 

1=1      ~U 
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Case  (a) : 

There  exists  a  feasible  solution  (h,Y)  having  h  =  h*  and 

1     ,   if  JL  =  e 

3%  ;   <  7*^    ,   if  &  =   m. 
0      ,   otherwise 

with  objective  function  value 

2  <  Z*  -  (C^  -  C.g) 

<  z   since  (c  .  -  c  „)  >  0. 
il     i8 

Hence,  there  is  a  contradiction  since  (h*,Y*)  cannot  be  an  optimal 

solution. 

Case  (b) : 

There  exists  a  feasible  solution  (h,Y)  having  h  =  h*  and 

1  ,   if  JL  =  6  -  [6/m.]m  and  I  >   0 

i   i 


y. .  +    [6/m. ]    ,   if  i  =   m. 
i  i 

,   otherwise 


with  objective  function  value 

z<z   -  (c.„  -  c     ) 

il     i,ra. 

l 

* 
<  z   since  (c.,  -  c    )  >  0  . 
il     i,m 
l 

Hence,  there  is  a  contradiction  since  (h*,Y*)  cannot  be  an  optimal 
solution. 

Q.E.D. 
The  optimal  algorithmic  solution  (h  ,Y°)  to  problem  P5.3  implies 
that  the  shortest  horizon  over  which  all  product  inventories  can  be 
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built  up  simultaneously  to  an  amount  I    is  h   production  periods. 

in 

ra. 

1 

During  this  horizon  each  product  should  be  assigned   E  y  .  times  where 

all  but  at  most  one  of  these  assignments  should  include  the  production 

m.-l 

1 

of  m.  batches  of  product  i.   If   E  y. .   equals  one,  then  the  other 

m.  -1 

i 

assignment  of  product  i  should  include  the  production  of   E  iy    . 

JU1      iX> 

batches.   Clearly,  since  problem  P5. 3  is  a  relaxed  version  of  problem 
P5.2,  (h  ,Y  )  is  also  an  optimal  solution  to  problem  P5.2  if  it  is 
feasible  to  problem  P5.2. 

Without  loss  of  generality,  it  can  be  assumed  that  the  first 

o 

y.     assignments  of  product  i  during  the  horizon  h  have  m  batches 
i,nti  i 

m. 

of  product  produced  where  E  y  .  assignments  are  made.   Let  z   denote 

JU1      U  ik 

the  number  of  these  machine-period  assignments  which  are  made  in  period  k. 

Then,  obviously,  (h  ,Y  )  can  be  converted  into  a  feasible  solution  to 

problem  P5.2  if  and  only  if  there  exists  a  feasible  solution  to  the 

following  constraint  set: 


h 

E     z        = 
k=l      lk 

l 

E     Yii 

k 

E     z        > 

<w..  /m.  ) 
ik      i 

N 

ik 

M 

1,2 N 


i  =  1,2,. . . ,N 
k  =  1,2 h°-l 

[P5.6] 

N 

1,2,.. . ,hv 
i=l 


'ik  S°'  integer  k=  l,2,...,hc 
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Since  E  y  ,  s  <w ,.°/m  >  for  all  i  =  1,2,...,N,  it  follows 
1=1      1X     xh        1 

directly  that  there  exists  a  feasible  solution  to  constraint  set  P5.6 
if  and  only  if 

N 

T,     <w.,  /m.  >  <  Mk 
.  -    lk  i 
i=l 

for  all  k  =  l,2,...,h  .   This,  coupled  with  the  following  result  and 

the  fact  that  any  feasible  solution  to  problem  P5.2  must  have  h  ^  h° 

(Property  2),  implies  that  problem  P5.2  has  a  feasible  solution  if  and 

only  if  (h  ,Y  )  can  be  converted  into  a  feasible  solution  to  problem  P5.2. 


Lemma  5.  3 


Problem  P5.2  has  a  feasible  solution  if  and  only  if 


E   ^wik/mi>  ==  ^  (5.14) 

i=l 


for  all  k  =  1,2,... ,h. 


Proof  of  Lemma  5.3 


The  "if"  part  follows  since  the  algorithmic  solution  (h°,Y°) 
has  already  been  shown  to  be  a  feasible  solution  when  (5.14)  is  satis- 
fied.  The  "only  if"  part  follows  by  summing  over  all  i  in  constraint 
(5.7),  assuming  m±   batches  are  produced  at  each  facility  assignment, 
and  summing  over  all  k  in  constraint  (5.8).   From  (5.7),  then,  any 
feasible  solution  must  satisfy 

xt    i     m-  B, 

N   k    i  i 

S        r        Z     Xi4i    *      Z   <Wii/V  <5-15> 

i=l   j=l    jtl  J         1=1      lk      x 

for   all   k  =    1,2,... ,h   and   from    (5.8)    any    feasible   solution  must    satisfy 
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N   k   m± 

Z   Z   E  xh  ;h  -  AIk  (5.16) 

i=l  j=l  jtl    J 

for  all  k  =  l,2,...,h.   Together,  relations  (5.15)  and  (5.16)  along 
with  the  integrality  requirement  on  x  .   imply  that  (5.14)  must  be 
satisfied  by  any  feasible  solution  to  problem  P5.2. 

Q.E.D. 
A  procedure  for  solving  problem  P5.2  is  to  solve  the  relaxed 
version  of  the  problem,  and  then  convert  this  solution  to  a  feasible 
solution  to  problem  P5.2  by  finding  a  feasible  solution  to  constraint 
set  P5.6.   Any  feasible  solution  determined  by  the  procedure  represents 
an  optimal  solution  to  problem  P5.  2  since  the  optimal  objective  function 
value  of  a  relaxed  version  of  a  problem  is  a  lower  bound  on  the  optimal 
objective  function  value  of  the  general  problem.   In  addition,  if  a 
feasible  solution  cannot  be  determined  to  constraint  set  P5.6,  then  it 
is  known  that  problem  P5.2  does  not  have  a  feasible  solution. 

Constraint  set  P5.6  describes  a  network  for  which  a  two-product, 
M-facility,  three-period  example  is  shown  in  Figure  5.2.   From  the  net- 
work it  can  be  seen  that  a  simple  procedure  for  determining  a  solution 

to  this  constraint  set  is  to  assign  each  product  i  (w   /m  )  times  to 

il   i 

vacant  facilities  in  period  one.   Then  assign  each  product  i  (w  n/m   )    - 

i2   1' 

<wil/mi>  times  to  vacant  facilities  in  periods  one  and  two,  and  in  gen- 
eral assign  each  product  i  <w  /m  >  -"<*.     /m)   additional  times  to 

IK    1  1 f K— J.    1 

vacant  facilities  in  periods  one  through  k  for  all  remaining  periods 
k  =  3,4,...,h  .   The  solution  obtained  is  optimal  unless  at  some  stage 
in  the  assignment  procedure  there  are  insufficient  facilities  available 
to  complete  all  assignments  in  which  case  problem  P5.2  does  not  have 
a  feasible  solution. 


Lower      so 
Bounds 


Upper 
Bounds 


Upper  and 
Lower  Bounds 

<Vmi> 


Figure  5.1   Network  Example  of  Constraint  Set  P5.6. 
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Computational  Experience 

The  algorithmic  procedure  was  programmed  in  Fortran  H  and  used 
to  solve  the  production  scheduling  problem.   An  IBM  370/165  computer 
was  used,  and  200  different  problems  were  considered  where  parameters 
were  generated  randomly.   Problems  having  as  many  as  20,000  variables 
and  20,000  constraints  were  generated.   All  problems  were  solved,  and 
the  computation  times  which  were  obtained  were  very  nearly  a  linear 
function  of  the  number  of  constraints.   The  average  solution  time  (as 
measured  in  terms  of  constraints  in  the  problem)  was  5.8  x  10 
seconds  per  constraint  with  a  range  of  approximately  ±  3  x  10~4 
seconds  per  constraint  about  the  mean.   The  computer  program  completed 
approximately  100,000  iterations  of  the  algorithm  each  second. 

Conclusions 

A  specially  structured  integer  program  has  been  identified, 
and  an  extremely  efficient  algorithm  for  determining  its  solution  has 
been  developed.   The  specially  structured  program  is  a  generalization 
of  a  multi-facility,  multi-product  production  scheduling  problem  which 
has  been  presented.   The  scheduling  problem  arises,  realistically,  in 
industrial  situations  where  schedulers  use  repetitive  cycle  plans  such 
as  the  well-known  economic  manufacturing  quantity  models  for  controlling 
production  systems. 


CHAPTER  6 
SUMMARY  AND  SUGGESTIONS  FOR  FUTURE  RESEARCH 

In  this  dissertation  a  class  of  multiple  facility,  multiple 
product,  production  inventory  scheduling  problems  has  been  considered. 
Specifically,  four  different  scenarios  for  this  class  have  been  addressed. 
In  a  straightforward  manner,  each  of  these  was  formulated  as  an  integer 
program.   It  was  then  shown  that  each  of  the  problems  could  be  reformu- 
lated as  a  network  flow  problem  which  can  be  solved  using  any  of  the 
very  efficient  algorithms  for  minimal-cost  flows  in  single  commodity 
networks.   In  addition,  for  some  important  special  cases,  single  pass 
procedures  were  developed  which  are  more  efficient  than  the  standard 
flow  algorithms.   In  the  fourth  scenario,  a  specially  structured  integer 
programming  formulation  was  identified,  and  an  extremely  efficient  algo- 
rithm was  developed  for  determining  its  solution. 

The  scenarios  which  have  been  considered  have  included  very 
general  and  realistic  assumptions.   The  originality  of  the  formulations 
for  this  class  of  problems,  in  addition  to  allowing  us  to  determine 
optimal  solutions  efficiently,  has  offered  considerable  insight  into 
the  operation  of  the  production  inventory  system.   These  insights  have 
allowed  us  to  enrich  the  problems  generously  to  include  many  other  real- 
istic aspects  of  the  production  system  heretofore  not  solvable  with 
present  techniques. 
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As  should  be  the  case,  however,  the  research  for  this  disserta- 
tion has  been  open-ended  in  character.   Each  question  answered  has 
allowed  new  questions  to  be  posed.   Some  of  the  logical  and  important 
extensions  which  it  would  seem  desirable  to  pursue  in  future  research  in 
this  area  include: 

1.  It  has  been  assumed  throughout  the  dissertation  that  setups 
occur  between  production  periods,  whereas  in  many  production  systems, 
setups  actually  occur  during  production  periods.   In  these  situations 
the  productive  capacity  of  the  system  is  depleted  when  setups  occur, 
and  an  important  extension  of  the  results  in  this  research  would  be  to 
obtain  equivalent  results  for  the  case  where  setups  are  assumed  to  occur 
during  the  production  period. 

2.  The  facility  usage  costs  which  have  been  used  do  not  take 
into  account  the  product  assigned  to  the  facility  during  the  previous 
period.  Realistically,  however,  the  usage  costs  often  should  only  be 
assessed  when  a  different  product  is  produced  on  the  facility  in  each 
of  the  two  periods.  Solving  the  production-inventory  scheduling  prob- 
lems for  this  more  general  case  would  be  a  considerable  contribution, 
but  the  problem  appears  to  be  quite  difficult. 

3.  Each  time  a  product  is  assigned  to  a  facility  for  one  period, 
it  has  been  assumed  that  only  an  integer  multiple  of  the  basic  batch 
size  for  each  product  can  be  produced.   An  important  generalization  of 
the  results  in  this  research  could  be  obtained  by  determining  equiva- 
lent results  for  the  case  where  partial  batches  are  allowed. 
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4.   Product  demand  rates,  though  not  necessarily  identical  in 
all  periods,  have  been  assumed  to  be  deterministic.   Considerable  addi- 
tional insight  into  the  multi-facility,  multi-product  production- 
inventory  system  could  be  achieved  by  studying  similar  problems  having 
stochastic  demand  rates. 


APPENDIX 


PROOF  OF  THEOREM  5.1 


Define 


h  =  TT  (in. p.) 
i=l   X  X 


and 


yU  = 


E  (d.,  /m.p.)     if   4  =  m. , 
lk   1  1  l 


otherwise 


i  =  1,2 N. 


(Al) 


Then,  from  (Al) , 


Jt=l    l£   k=l   lk   * 


=   w.,  -  (I   -  I.n)/P. 
lh     iH    iO    x 


for  all   i  =  1,2,. . . ,N,  and 


N   mi        h    N 


E   E  y  .  =   E   E  (d  /m  p  ) 
i=l  A=l   lJi   k=l  i=l   lk  X  1 


<   E  M,   since  from  (5.11) 
k=l 


L    (d   /m.p. )  <  M. 
.  ,       lk   i  i 
i=l 


(A2) 


Relationship  (A2)  implies  that 


Mh  -  E   E  y  .  =  s  s  1. 
i=l  X=l   1 
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Thus,  during  the  horizon  h,  all  product  demands  are  satisfied 
exactly,  and  there  are  s  additional  unassigned  periods  available. 
Hence,  by  assigning  each  product  i,  (j,.  additional  times  during  the 
horizon  h  where 


0  <  u.  <  max  [0  ,  (T^  -  I1())/P1>  } 


and 


N 

E  n  =  s, 

i=l  X 

at  least  one  product  must  become  discretely  closer  to  satisfying  con- 
straint (5.9).  If  all  constraints  are  now  satisfied,  then  a  feasible 
solution  has  been  determined.   Otherwise,  redefine 

(fiH  -  fiO)/Pi  =  (fiH  "  TiO)/Pi  -  ^iPimi 
and 

h+1,  h+2,  ...  =  1,2 

and  repeat  the  procedure.   A  feasible  solution  must  be  obtained  in 
a  finite  number  of  iterations  since  at  each  iteration,  all  demands 
are  satisfied  and  at  least  one  product  is  assigned  at  least  one 
additional  time. 

Q.E.D. 
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